N PS  ARCHIVE 

1958 
FAHLAND,  F. 


RvlkHi 


ANALOGY 


S  OF  FLOW  THROUGH  CASCADES 


FRANK  R,  FAHLAND 


'It1 


t\*M 


■  ■  ■ 

Library 

Naval  Postgraduate  Sohool 

Monterey,  California 


AN  ELECTRICAL  ANALOGY 
FOR 
ANALYSIS  OF  FLO'"'  THROUGH  CASCADES 

*  *  *  * 

Frank   R.    Fahland 

and 
Larry   L.    Havkins 


AN  ELECTRICAL  ANALOGY 
FOR 
ANALYSIS  OF  FLOV  THROUGH  CASCADES 


.  by 

Frank  R.  Fahland 
il 

Lieutenant,  United  States  Navy 
and 
Larry  L.  Hawkins 
Lieutenant,  United  States  Navy 


Submitted  in  partial  fulfillment  of 
the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 
IN 
MECHANICAL  ENGINEERING 


United  States  Naval  Postgraduate  School 
Monterey,  California 

19  5  8 


F5   fWC 


Library 

Naval  Postgraduate  School 

Monterey,  California 


AN  ELECTRICAL  ANALOGY 
FOR 
ANALYSIS  OF  FLOW  THROUGH  CASCADES 

by 

Frank  R.  Fahland 

and 
Larry  L.  Hawkins 


This  work  is  accepted  as  fulfilling 
the  thesis  requirements  for  the  degree  of 
MASTER  OF  SCIENCE 
IN 
MECHANICAL  ENGINEERING 

from  the 
United  States  Naval  Postgraduate  School 


SUMMARY 

A  method  of  analysis  of  arbitrary  cascades  under  conditions 
of  ideal  fluid  flow  is  described.   The  method  is  based  on  con- 
formal  transformation,  achieved  experimentally  through  an  elec- 
trical analogy.   Tests  were  conducted  to  permit  comparison  of 
results  obtained  by  this  method  with  those  of  a  previous  ana- 
lytical solution  and  with  NACA  two  dimensional  cascade  tests. 
Results  indicate  that  the  quick  and  inexpensive  procedure  of 
this  method  is  applicable  to  engineering  problems  of  cascade 
analysis . 

This  work  was  conducted  at  the  U.  S.  Naval  Postgraduate 
School,  Monterey,  California. 
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1.   INTRODUCTION 

One  of  the  first  steps  in  the  aerodynamic  design  of  an 
axial  flow  turbomachine  is  the  design  of  cascades  which  produce 
desired  deflections  of  the  fluid  as  it  passes  through  the  ma- 
chine.  Fundamental  to  the  solution  of  this  problem  is  the  de- 
termination of  the  two  dimensional  flow  of  an  ideal  fluid  across 
cascades  of  compressor  and  turbine  blades.   An  ideal  fluid  is 
here  assumed  to  mean  one  that  is  homogeneous,  continuous,  incom- 
pressible, and  inviscid. 

The  development  of  a  method  which  predicts  the  two  dimen-. 
sional  flow  of  an  ideal  fluid  across  any  arbitrary  cascade  was 
the  objective  of  this  investigation.   One  measure  of  cascade 
performance  is  the  fluid  turning  angle.   More  basic  parameters 
are  the  cascade  coefficient  K    defined  in  Section  5.1  and  the 
relative  zero  lift  angle  6    defined  in  Section  5.1  and  illus- 
trated  in  Fig.  17.   Of  primary  interest  was  the  effect  of  changes 
in  cascade  solidity  and  blade  stagger  on  these  indicators  of  cas- 
cade performance,  but  a  preliminary  test  of  the  ability  of  this 
method  to  predict  velocity  distributiDn  around  the  airfoils  in 
cascade  is  also  included. 

Considerable  ingenuity  has  been  exerted  in  previous  work 
done  on  this  problem,  and  many  mathematical  analyses  exist  of 
the  two  dimensional  potential  flow  past  infinite  cascades  of 
airfoils.   Some  provide  exact  solutions  for  special  shapes  such 
as  straight  blades  of  zero  thickness.   See  Ref s .  1  and  2.   Some 
provide  approximate  solutions  for  particular  airfoil  shapes  and 


cascade  arrangements.   See  Refs.  3  and  4.   Finally  in  regard  to 
the  general  problem  of  arbitrary  airfoils  in  cascade,  there 
have  been  several  attempts  at  a  mathematical  solution.   These 
usually  involve  iteration  methods.   See  Refs.  5,  6,  7,  and  8. 
The  labor  involved  in  the  application  of  these  methods  is  great; 
and,  for  some  applications,  even  modern  high  speed  computers 
prove  inadequate.   Particular  difficulty  is  encountered  with 
blades  of  high  camber  and  cascades  of  high  stagger  angles.   Hence 
the  engineering  application  of  these  methods  is  limited. 

Theoretically,  the  method  here  investigated  is  based  upon 
the  conformal  transformation  of  an  arbitrary  cascade  into  a  cir- 
cle.  This  is  accomplished  experimentally,  using  electrical  con- 
ducting paper.   The  flow  of  an  ideal  fluid  is  analytically  known 
about  a  circle.   Through  Laplace's  equation  the  flow  can  be  de- 
termined about  the  cascade. 

To  compare  the  results  obtained  by  this  analog  method  with 
those  of  a  previous  analytical  study,  the  airfoil  in  cascade 
which  had  been  analyzed  by  Mutterperl  and  reported  in  Ref .  7 
was  investigated.   To  determine  the  effect  of  solidity  and  stag- 
ger upon  turning  angle,  a  blade  more  closely  approximating  those 
currectly  used  in  the  gas  turbine  industry  was  desired.   The 
NACA  65-(15)10  airfoil  section  was  chosen  for  this  reason  and 
because  its  use  permitted  comparison  with  NACA  two  dimensional 
cascade  tests  reported  in  Ref.  9. 

This  work  was  conducted  at  the  U.  S.  Naval  Postgraduate 
School,  Monterey,  California  during  the  1957  -  1958  school  year. 

The  authors  wish  to  express  their  appreciation  to  Professor 
M.  H.  Vavra  and  Professor  T.  H.  Gawain  of  the  U.  S.  Naval  Post- 


graduate  School.   This  thesis  project  was  conceived  by  Professor 
Vavra  and  is  an  outgrowth  of  previous  investigative  work  per- 
formed by  Professors  Vavra  and  Gawain.   The  association  with 
Professor  Gawain  who  guided  and  assisted  in  the  work  of  this  in- 
vestigation was  most  rewarding. 


2.   TABLE  OF  SYMBOLS 

a       radius  of  circle  in  the  picture  plane,  feet 

b  distance  from  origin  of  sources  and  sinks  in  picture 

plane,  feet 

/C       chord  length,  inches  or  feet 

/Ct  airfoil  section  coefficient  of  lift 

I  dimensionless  transformation  parameter,  defined  by 

equation  3.27 

I  current,  amperes 

I  indicates  imaginary  part  in  complex  number  notation 

H  characteristic  distance  used  in  conformal  transforma- 

tion, feet 

m       source  strength  in  picture  plane,  feet  squared  per  sec- 
ond 

q  dynamic  head,  pounds  per  square  foot 

r  radius  vector  length  in  picture  plane,  feet 

s  tangential  spacing  betv/een  blades,  inches  or  feet 

w  complex  velocity  potential  function 

x  real  axis,  real  plane 

y  imaginary  axis,  real  plane 

z  complex  coordinate,  real  plane 

E  potential,  volts 

I  total  current  function,  defined  in  Appendix  A 


K  cascade  coefficient,  defined  in  Section  5.1 

R  total  resistance,  ohms   /  ^ 

S  pressure  coefficient  =  (  — 

V  velocity,  feet  per  second 


in 


oC  angle  of  fluid  flow,  degrees 

oCQ         absolute  angle  of  attack 

ofm   angle  of  mean  velocity 

of.    angle  of  zero  lift 

oCl         cascade  entrance  angle 

oCz         cascade  exit  angle 

oQ    ideal  approach  angle 
P  cascade  stagger  angle,  degrees 

6L       relative  angle  of  zero  lift,  degrees 

r  ioge  £- 

$  complex  coordinate  in  the  picture  plane 

n  imaginary  axis,  picture  plane 

Q  angular  coordinate,  picture  plane,  degrees 

V  specific  resistance,  ohms 

|  real  axis,  picture  plane 

p  resistivity,  ohm-inches 

p  density,  slugs  per  cubic  foot 

(T  solidity  =  -y- 

X  'aspect  ratio',  defined  by  equation  3.13 

f  velocity  potential  function 

v,  velocity  stream  function 

I  circulation,  feet  squared  per  second 

Q  turning  angle,  degrees 

A  vortex  strength,  feet  squared  per  second 

u  percent  potential 

^  percent  current  quantity 


Subscripts : 

1  upstream  of  blade  toy.' 

2  dovnstream  of  blade  row 

4 

a  axial  direction 

i  ideal 

m  mean 

s  stagnation 

t  tangential  direction 

L  leading  edge 

T  trailing  edge 
Superscripts : 

*  indicates  maximum  value 


3.   ANALYSIS 
3.1   IDEAL  FLUID  FLO".' 

If,  in  a  two  dimensional  steady  flow,  there  exists  a  veloc- 
ity point  function  such  that 


V  =  v(*)=V(*>i/) 


3.1 


the  necessary  and  sufficient  condition  for  this  to  be  potential 
flov  is  that 

VX\7  =  0  3.2 

which  means  that 

y^vf  3.3 

But,  if  further,  the  density  of  this  flo1"  be  constant,  continu- 
ity requires  that 

V*V  =  3  3.4 

If  both  3.3  and  3.4  are  satisfied;  then  the  potential  is  a  plane 

0 

harmonic  function  such  that 

V>=0  3.5 

which  is  LaPlace's  equation. 

Correspondingly  in  this  theoretical,  two  dimensional  fluid 
flow  there  is  a  stream  function  for  which 

V>=0  3.6 


Trajectories  of  constant  <p     are  orthogonal  to  trajectories  of 
constant  ~\p  . 

Analagous  to  this  theoretical,  two  dimensional  fluid  flow 
is  the  flow  of  a  steady  current  in  a  thin,  flat  sheet  of  uni- 
form conductivity.   Details  of  the  analogy  are  contained  in  Ap- 
pendix A.   On  the  sheet  the  electric  potential  E,  analagous  to 
(D   ,  can  be  represented  by 

Similarly,  introduction  of  a  total  current  function  1  results  in 

VV  =  0  3.8 

Trajectories  of  constant  E  are  orthogonal  to  trajectories  of 
constant  1 . 

3.2  CONFORMAL  TRANSFORMATION 

Various  flows  of  an  ideal  fluid  are  analytically  known  a- 
bout  certain  simple  domains  such  as  a  circle.   Laplace's  equa- 
tion, which  characterizes  such  a  perfect  fluid  flow,  remains  in- 
variant for  any  conformal  transformation.   Hence,  if  an  arbitrary 
shape  can  be  conformally  transformed  into  a  suitable  circle,  the 
flow  of  an  ideal  fluid  about  the  arbitrary  shape  can  be  deter- 
mined . 

The  theory  of  conformal  transformation  with  the  aid  of  an 
electric  tank  is  discussed  in  Ref .  10.   In  general  the  whole  re- 
gion between  two  arbitrarily  specified  closed  curves  is  not  trans- 
formable into  the  whole  region  between  two  other  arbitrarily 
specified  closed  curves.   Transformation  can  be  accomplished, 
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however,  if  there  exists  a  function,  in  the  picture  plane,  of 
the  real  plane  complex  variable  2  such  that  for  every  assigned 
value  of  the  variable  the  function  has  a  definite  value  or  sys- 
tem of  values  and  also  a  definite  derivative. 

The  arbitrary  shape  transformed  in  this  investigation  con- 
sists of  an  infinite  cascade  of  equally  spaced  compressor 
blades  of  the  same  shape  and  orientation. 

Imagine  a  conducting  sheet  of  uniform  conductivity,  infin- 
ite in  extent  in  two  dimensions,  Fig.  1,  in  which  the  blades 
are  areas  of  infinite  resistance  (cutouts  from  the  sheet).   With 
a  potential  E*  imposed  on  a  line  parallel  to  the  axis  of  the 
cascade  and  located  infinitely  far  ahead  of  the  cascade  and  a 
zero  potential  imposed  on  a  similar  line  located  infinitely  far 
behind  the  cascade,  corresponding  points  on  each  of  the  blades 
have  equal  potential.   Immediately  in  front  of  the  cascade  a 
line  connecting  points  of  equal  potential  exhibits  a  shape  that 
is  strongly  influenced  by  the  geometry  of  the  cascade;  but  with 
increasing  separation  from  the  cascade  the  influence  diminishes 
rapidly,  and  at  finite  distances  ahead  of  and  behind  the  cascade 
these  lines  of  constant  potential  are  virtually  straight  and 
parallel  to  the  axis  of  the  cascade.   Hence,  the  conditions  at 
infinity  may  be  closely  approximated  by  imposing  a  voltage  on 
the  sheet  through  straight  bus  bars  parallel  to  the  axis  of  the 
cascade  and  a  finite  distance  away. 

No  generality  is  lost  if  a  single  strip  of  width  s  which  in- 
cludes the  space  between  two  blade  cutouts  is  considered,  provid- 
ed that  the  boundaries  or  edges  of  the  strip  are  lines  of  con- 
stant 1 .   Let  the  boundaries  be  chosen  to  coincide  with  portions 


of  the  blade  profiles  as  shown  in  Fig.  2. 
For  the  strip  of  conducting  sheet 

R*i*=E*  3.9 

For  points  on  this  strip  (from  Appendix  A) 

E+ivI   =f(2)  3.10 

Dividing  both  sides  of  3.10  by  3.9 

{,(«)*£,+!££-,  =6tU  3.ii 

\*  will  be  constant  along  the  edge  of  the  strip  and  £  will  vary 
from  one  on  the  front  bus  bar  to  zero  on  the  after  bus  bar.  Any 
point  on  the  strip  can  be  located  by  coordinates  x,  y;  but  O  and 
Vdir  are  continuous  functions  on  the  str,ip  and  are  definite  at  any 
point.  Optionally,  then,  any  point  may  be  designated  by  its  co- 
ordinates or  by  the  values  of  o  and  <d*  which  characterize  it. 
In  symbols : 

and  conversely  3.12 

On  the  front  part  of  the  profile  cut-out  G  will  reach  a  maxi- 
mum value  at  a  point   1,  and  on  the  trailing  part  a  minimum 
value  at  a  point  2.   Fig.  3  shows  a  typical  plot  of  percent  po- 
tential £/  vs.  percent  chord.   The  quantity  Q,   —  o:  is  one  of 
the  parameters  used  to  accomplish  transformation. 

In  the  region  near  the  bus  bars  where  lines  of  constant  o 
are  straight  and  parallel  to  the  cascade,  imagine  two  lines  of 

10 


constants  ;  Ga  and  £.      such  that  Ga— £/a,  =  U(  ""  Sz   '      See  Fi9«  4« 
Imagine  that  the  width  of  the  conducting  strip,  s,  is  divided 
into  T)t    equal  parts,  each  measuring  -~-  .   This  same  increment 
of  length  is  used  to  measure  along  the  distance  -i.  between  C/a 

andCiv»  dividing  this  distance  into  7)i=  -£- parts.   Lines  of 

constant  £  drawn  at  these  stations  are  separated  by  equal  dis- 
tances and  represent  equal  increments A£,  .   For  convenience  it 
is  assumed  in  Fig.  4  that  f),  and  7\  are  integers.   Lines  of 
constant  G  and  constant  w  drawn  at  these  measured  stations 
form  an  orthogonal  grid  of  T\f\x      squares  over  this  portion  of 
the  strip. 

A  similar  conceptual  procedure  is  valid  at  any  position  in 
the  strip.   Lines  of  constant  Q,     may  be  drawn  such  that 
£■  ""*Cj  -  U,  —  C/^  .   Let  the  distance  between  them,  as  measured 
along  a  line  of  constant  1,  be  divided  into  T)z  segments  of  equal 
potential  drop,  Aq  .   The  distance  between  the  edges  of  the 
strip,  as  measured  along  a  line  of  constant  £  »  can  then  be  di- 
vided into  T),  segments.   A  grid  of  orthogonal  curvilinear  co- 
ordinates results  and  the  region  consists  of  the  same  number 
T^TYj  of  curvilinear  squares  as  before. 

Let  there  be  a  parameter  T  such  that 

In  the  regions  of  the  bus  bars  G  varies  linearly  with  1;  so 

(-Vf)„ s  6«JCii  3-14 

0 

V    then  reduces  to  the  ratio  ■»-  which  might  be  considered  an 

o 

'aspect  ratio'  for  the  gridded  section.   It  represents  the  ra- 
il 


tio  of  the  number  of  squares  along  the  strip  to  the  number  of 
squares  across  the  strip.   The  quantity  represents  a  unique  and 
invariant  characteristic  of  the  geometry,  and  must  be  maintained 
in  mapping  to  achieve  conformal  transformation. 

Consider  the  denominator  of  the  quantity  V  .   From  the  anal- 
ogy developed  in  Appendix  A,  at  any  point, 

In  the  region  near  the  bus  bars  where  \        is  constant  across  the 
width  of  the  strip 

I 

Multiplying  both  sides  by  s,  the  width  of  the  strip, 

Vl*S    =  (Ea-Eb)(f)'  3.17 

• 

<?2"5  is  the  total  current  flowing  in  the  strip  between  bus  bars. 
The  right  side  of  the  equation  might  be  labeled  the  specific 
voltage  drop  of  the  strip.   It  represents  the  voltage  drop  be- 
tween two  equipotential  lines  which  form  a  square  with  lines  of 
constant  I   that  are  s  distance  apart.   Equations  3.15,  3.16,  and 
3.17  may  be  used  in  non-dimensional  form  if  each  is  divided  by 
equation  3.9. 

Now  consider  the  %     or  picture  plane  of  Fig.  4  with  sources 

a* 
of  equal  strength  m  on  the  negative,  real  axis  at  -b  and  — — 

fa 

and  sinks  of  the  same  strength  on  the  positive,  real  axis  at  +b 

and  +  -=-   .   These  sources  and  sinks  create  a  complex  function 
b 


in  the  fc  plane.   See  Appendix  B: 


ur  =  -**- 

ZTT 


-     <p+  If 


3.18 


Any  point  in  the  picture  plane  can  be  located  by  coordinates  T], 
£  or  by  the  values  of  p   and  ~W  which  characterize  the  point. 
Symbolically, 


>=y(^7l) 


and  conversely  3.19 


The  circle   5   :CL  is  a  line  of  constant  1M   .   Sought  is  the 
conformal  transformation  of  the  strip  of  the  real  plane  external 
to  the  blade  cut-out  onto  the  area  of  the  picture  plane  outside 
the  circle  so  that  the  cascade  blade  outline  becomes  the  circum- 
ference of  the  circle.  O   is  a  minimum  at  a  point  1  and  a  maxi- 
mum at  a  point  2.   The  difference  in  potential  p  —  <D         is 


3.20 


TT        cosh  T- 1  m 


3.21 


See  Appendix  B.   It  will  be  shown  that  the  right  side  of  eq.  3.21 
is  the  corresponding  form  of  f     for  the  picture  plane. 

With  reference  to  the  strength  term  m  in  the  denominator  of 
eq .  3.21,  recall  that  in  two  dimensional,  unidirectional  flow  of 
an  ideal  fluid  the  uniform  quantity  rate  of  flow  between  stream- 
lines h  distance  apart  is  Vh,  Fig.  5.  Correspondingly  for  a 
point  source  from  which  flow  emanates,  a  uniform  quantity  rate 
of  flow  crosses  any  simply  connected  curve  enclosing  the  source. 
The  velocity  at  any  point  is  the  gradient  of  p    ,    that  is  V  =  yy»~» 
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The  change  of  p      between  two  <P     countours  which,  with  two 
streamlines  h  distance  apart,  form  a  curvilinear  square  is 
— L—h   .   This  quantity  might  be  dubbed  the  specific  p      drop, 
analogous  to  specific  voltage  drop,  but  for  fluid  flow  consid- 
erations it  is  designated  a  quantity  rate  of  flow,  here  repre- 
sented by  m. 

To  achieve  conformal  transformation  between  the  z  and  5* 
planes,  the  following  condition  is  imposed 

P^K,-kz&  3.22 

where  /<,  and   /<z  are  arbitrary  constants. 

It  follows  that  <Pz~f,    =  kz  (£,-£*.)   .  3.23 

Further,  let   TT)  -   /<25£-~-J   ,  3.24 

Then       %"fl  sJbJ&LdLL  3.25 

The  consequence  of  equal  V     in  the  real  and  picture  planes  is 
this:   for  any  area  in  the  real  plane  bounded  by  lines  of  con- 
stant <Z/    and  w  there  is  an  area  in  the  picture  plane  bounded  by 
corresponding  lines  of  constant  p     and  "\js    .   The  'aspect  ratios' 
of  the  areas  are  the  same.  See  Fig.  4. 

For  every  point  o  ,  ^ir  in  the  real  plane  there  is  a  cor- 
responding point  y  ,  Y'   in  the  picture  plane.  — —    can  be  meas- 
ured in  the  real  plane,  and  a  corresponding  -=£-    can  be  calcu- 
lated  in  the  picture  plane.   These  conditions  are  sufficient  to 
establish  the  existence  of  a  function  of  a  complex  variable  which 
relates  the  two  planes  conformally.   An  analytical  expression  of 
this  function  is  not  determined  nor  required. 
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Definition  of  V     in  the  5   plane  uniquely  determines 
the  location  of  the  sources  and  sinks.   Eq .  3.21  is  solved  giv- 
ing 

irt 

cosh  If   =  J? ±_L  o  9A 

Further,  the  periodicity  of  the  mapping  function  is  such  that 
each  strip  of  width  s  in  the  infinite  cascade  is  mapped  onto 
this  same  area  in  the  picture  plane;  so  the  result  is  truly  a 
map  of  the  infinite  cascade.   Fig.  6  indicates  the  location  of 
corresponding  points  on  the  two  planes. 

Reference  to  eqs.  3.23,  3.24  and  3.25  indicates  that  only 
three  quantities  need  oe  measured  in  the  real  plane  to  accom- 
plish the  transformation.   These  are  fc,""£»z  j  \^x)^o  i       5. 
Other  measurable  quantities  might  perhaps  be  used,  but  the  a- 
bove  listed  parameters  seem  to  be  the  ones  most  naturally,  con- 
veniently, and  accurately  obtainable.   They  are  sufficient  to 
achieve  transformation. 

The  position  of  a  point  on  the  blade  profile  such  as  P  in 
Fig.  6  may  be  expressed  in  terms  of  the  measured  value  of  C 
which  characterizes  it.   Similarly  the  position  of  a  point  on 
the  circle  such  as  P'  in  Fig.  6  may  be  expressed  in  terms  of  the 
calculated  value  of  <P  which  characterizes  it.   Points  P  on  the 
blade  and  P'  on  the  circle  which  correspond  have  identical  val- 
ues of  the  parameter  f,  where  f  is  defined  by  tie  expression 

Each  such  point  P'  on  the  circle  also  has  an  angular  co- 
ordinate Q   as  shown  in  Fig.  7.   The  coordinate  Q   is  related  to 
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the  parameter  f  by  the  expression 


Q   =  cos" 


cosh  7f   -£- 


-  / 


c*t-(*-o    , 


3.28 


This  formula  is  derived  in  Appendix  B. 


3.3   THE  FLOW  RELATIONS 


Once  the  transformation  has  been  completed  the  flow  in  the 
picture  plane  can  be  modified  to  correspond  to  various  desired 
flow  conditions  in  the  real  plane.   This  involves  variation  in 
the  direction  of  the  mean  flow  and  variation  in  the  circulation 
around  the  blade  to  satisfy  the  Kutta  condition. 

The  direction  of  the  mean  flow  may  be  changed  by  super- 
position of  a  set  of  vortices  j\   at  the  singularities  in  the 
picture  plane  as  shown  in  Fig.  7.   The  circulation  around  the 
airfoil  can  be  modified  by  superposition  of  a  second  set  of 
vortices  -I—  as  indicated  in  the  same  figure. 

The  Kutta  condition  states  that  the  rear  stagnation  point 
of  the  fluid  flow  about  an  airfoil  with  a  sharp  trailing  edge 
occurs  at  the  sharp  trailing. edge.   Therefore  the  rear  stagna- 
tion point  in  the  picture  plane  must  occur  at  the  point  on  the 
circle  which  corresponds  to  the  trailing  edge.   This  determines 
the  strength  of  the  circulation  that  must  be  imposed. 

At  the  desired  angle  of  attack  and  with  the  Kutta  condition 
satisfied,  the  velocity   v$  at  every  point  on  the  circle  may  be 
calculated.   Since  the  conformal  transformation  has  established 
the  correspondence  between  points  on  the  blade  in  the  real  plane 
and  on  the  circle  in  the  picture  plane,  the  velocity  V2  at  each 
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point  on  the  blade  is  obtainable  by  means  of  the  relation 

It  may  further  be  shown  that4r  =  4l — =  fi   1  ffi~"  ffi  ) 

dS   df     _d£   \5/-<S,  / 

d2      dy 


3.29 
3.30 


where 


A< 


A2 


is  measured  along  the  blade  profile.   See  Appendix  B. 


It  should  be  noted  that  when  transformation  has  been  accomplished 
velocities  of  corresponding  points  on  the  blade  and  circle  are 


related  by  the  factor 


di 


d* 


,  which  is  a  constant  at  any  given 
point  and  is  not  a  function  of  angle  of  attack  or  circulation  im- 
posed upon  the  circle  in  the  picture  plane. 

The  detailed  analytical  development  of  the  ideal  fluid  flow 
in  the  picture  plane  is  contained  in  Appendix  C.   What  follows 
here  is  a  descriptive  summary  of  the  procedure.   The  sources  m 
placed  in  the  picture  plane  to  accomplish  the  conformal  trans- 
formation create  a  'straight  through'  fluid  flow.   This  flow 
causes  the  front  and  rear  stagnation  points  to  occur  on  the  cir- 
cle at  the  ends  of  a  diameter  which  lies  along  the  real  axis. 
See  points  1  and  2  in  Fig.  7. 

To  develop  an  angle  of  attack  of  the  fluid  with  the  zero 
direction  of  the  real  axis,  vortices  of  strength -A.  are  imposed 
upon  the  four  point  sources.   For  positive  approach  angles,  coun- 
ter-clockwise vortices  of  equal  strength  are  placed  at  the  source 

at  £  =  —  b  and  at  the  sink  at  £  —  ~r-       •   Clockwise  vortices  of 

b 

the  same  strength  are  placed  at  the  other  source  and  sink.   The 
mean  approach  angle  then  becomes: 


0C™  =  tan 


-i  A 


m 


3.31 
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The  front  and  rear  stagnation  points  still  occur  at  the  ends  of 
a  diameter  of  the  circle  but  the  diameter  makes  an  angle  with 
the  real  axis 

6s  =  tern"  (tan  </mt<mhr)  3#32 

The  flow  is  undeflected  as  it  passes  around  the  circle;  so  no 

net  circulation  or  lift  is  present. 

To  produce  deflection  and  lift,  vortices  each  of  strength 

-4-  and  of  the  same  direction  are  placed  at  the  source  and  sink 

inside  the  circle.   The  strength  and  the  direction  of  these 

vortices  are  those  required  to  satisfy  the  Kutta  condition; 

i.e.  the  rear  stagnation  point  must  be  brought  to  the  location 

on  the  circle  which  corresponds  through  mapping  to  the  trailing 

edge  of  the  airfoil.   Vortices  of  the  opposite  direction  must 

be  placed  at  the  source  and  sink  outside  the  circle  to  insure 

that  the  stagnation  streamline  remains  a  circle  of  radius  0. . 

With  sources  and  sinks  m.  and  with  vortices  -A-  and  -I— 

Z 

in   place,    the   picture   plane   is   complete.      The  velocity   at   any 
point   on   the   circle   is 


LA    _   w  „  -zm  (coshrsinQ      \+2A.(stnhXiosQ  __\- T  (  sinhzi  \ 

V$LrLa  VS ~     TTQ.  \co5h2r-5?7U<9/    rra[coshlZSinZ&J  2tra\coshZ1-&riZG  )  3  ^ 
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m  is   determined   by     777  =  s(Va)  =  Sv^r"/^  3.34 

-~-     is   determined   by   eq .    3.31   once    aCm    is   arbitrarily  chosen. 
To   satisfy   the  Kutta    condition,    let  8=  QT   and        vq  =0.       The  val- 
ue  of    Qj.    is   determined   by   substituting    the   value   of      fj,         in 
eq.    3.28.      Solving   eq.    3.33   then  gives 

F      /^y      cos  &r        sin  6r 
-£-=  lanxm-crtFt  ~  sinhr  3.35 
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The  picture  plane  is  thus  completely  determined.   See  Fig.  7. 
For  subsequent  work,  formulas  for  particular  quantities 
relating  to  the  flow  direction  are  developed.   The  cascade  en- 
trance angle  is 


aC,  =  tan 


-i 


tarioCm  +  J_ 
Zm 


3.36 


The  cascade  leaving  angle  is 


aC  -  Tan' 


-> 


tan  oTm  -  V 


Zm 


3.37 


The  turning  angle  of  the  fluid  as  it  passes  through  the  cascade 
is 


6=°c,-oCz 


3.38 


Of  aerodynamic  interest  is  the  ideal  angle  of  attack.   When  the 
front  stagnation  point  occurs  at  the  point  of  the  circle  corres- 
ponding to  the  leading  edge  of  the  airfoil 


o<i  =  o(,  =  tan"    s»QT-sir>ej, 

_tQnh^9r-co5et)_ 


3.39 


The  lift  on  any  blade  is  developed  from  the  elementary  relation 


L=fV™r  =  *,qA 


pc     lX 


3.40 


p   is  the  density  of  the  fluid  and  f\    is  the  projected  area  of 
the  section  along  the  chord.   The  subscript  X  appearing  in  the 
above  listed  equations  indicates  that  the  coefficient  of  lift 
and  the  dynamic  head  c|  may  be  based  upon  the  cascade  entrance 
velocity  \J{     ,    the  cascade  leaving  velocity  \/i     ,  or  the  cascade 
mean  velocity  Y7T)  •   Each  of  the  three  coefficients  of  lift  has 
its  reasons  for  being  used  in  the  analysis  of  cascade  performance. 
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/Co   ,    based  upon  Vm  ,  permits  an  aerodynamic  tie-in  with  single 
airfoil  theory,  /Co      ,    based  upon  V,      is  useful  in  cascade  design 
because  conditions  upstream  of  the  cascade  are  generally  known. 
Co      based  upon  V2   ,  makes  allowance  for  the  deceleration  of 
the  flow  as  it  passes  through  a  compressor  cascade.   Ref.  11 
states  that  acutal  cascade  tests  have  shown  that  the  drag-lift 
ratio  reaches  a  minimum  at  an  approximately  constant  value  of 

In  terms  of  the  quantities  determined  by  the  previous  e- 
quations,  the  coefficient  of  lift  is 

^*-£sffe(ta»*  -*»*"*)  3.41 
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4 .   EQUIPMENT 

The  equipment  consisted  of  the  specially  designed  apparatus 
used  to  hold  the  two  dimensional,  conducting  paper,  with  appurte- 
nances to  impose  the  potential  and  complete  the  circuit.   Associ- 
ated with  this  piece  of  equipment  were  various  electrical  in- 
struments necessary  for  potential  supply  and  measurements.   A 
list  of  the  instruments  with  specifications,  where  applicable, 
is  given  in  Table  1 . 

A  complete  test  setup  is  shown  in  Fig.  8.   The  figure  shows 
the  conducting  paper,  the  working  surface  and  the  associated  e- 
lectrical  instruments.   The  remaining  item  is  the  profile  tem- 
plate used  as  a  guide  in  cutting  the  conducting  paper. 

4.1   THE  WORKING  SURFACE 

Fig.  9  shows  the  original  working  surface  which  existed  at 
the  start  of  the  investigation.   It  was  constructed  of  3/4  inch 
plywood  overlaid  with  3/16  inch  Masonite  (rough  surface  toward 
the  plywood)  and  was  fitted  with  Birmingham  Gage  Number  25  copper 
bus  bars,  one  inch  wide  and  five  feet  long.   These  bus  bars  were 
positioned  on  the  table  parallel  and  29.5  inches  apart  at  their 
inside  edges.   They  served  as  the  electrodes  for  imposing  a  po- 
tential difference  across  the  paper  cascade. 

The  paper  was  placed  over  the  bus  bars  with  its  conducting 
(graphite)  surface  toward  the  Masonite.   It  was  clamped  as 
shown  in  Fig.  10.   This  method  of  clamping  provided  a  tight, 
line  contact  between  the  bus  bars  and  the  paper. 
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This  setup  left  the  ends  of  the  paper  as  unknown  boundaries 
which  were  streamlines  of  the  flow.   It  was  necessary  to  use  an 
iteration  process  to  determine  these  boundaries.   A  first  approx- 
imation was  made  by  eye  and  the  paper  cut  accordingly.   The  suc- 
ceeding approximations  were  then  obtained  by  flux  plotting  a 
streamline  from  experimentally  determined  constant  potential 
lines.   These  plots  were  made  near  the  middle  of  the  cascade  to 
minimize  effects  of  the  previous  boundaries.   The  cascade  ends 
were  then  cut  to  this  nev ly  determined  streamline.   This  process 
of  refinement  was  continued  until  the  correction  required  became 
as  small  as  the  errors  of  measurement. 

In  order  to  by-pass  the  time  and  labor  consumed  in  obtaining 
these  boundaries  a  second  table  was  constructed.   This  table 
shown  in  Fig.  11,  was  similar  to  the  original  in  its  method  of 
connecting  electrodes  to  the  conducting  paper.   The  contact  be- 
tween the  electrodes  and  the  paper  was  obtained  by  the  same  man- 
ner of  clamping.   The  copper  bus  bars  were  widened  to  32  inches 
however,  and  maintained  parallel  to  +  0.003  inches  in  a  three 
foot  length. 

The  table  was  lengthened  to  six  feet  overall  with  a  foot 
and  a  half  of  each  end  hinged  to  the  main,  three  foot  section, 
"ith  this  provision,  the  paper  was  clamped  to  the  table  as  in 
Figs.  10  and  12  and  the  ends  were  folded  down  under  the  main 
section  thereby  bringing  the  ends  of  the  paper  cascade  together 
in  a  physical  and  therefore  electrical  contact.   The  cascade 
then  had  no  beginning  or  end  and  was,  in  effect,  infinite. 

Since  cascades  of  varying  solidity  required  various  lengths 
of  paper  in  order  to  come  to  contact  at  the  joint,  it  vas  neces- 
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sary  only  to  lengthen  or  shorten  the  bight  of  the  paper  at  the 
hinges  when  securing  the  paper  to  the  table  prior  to  folding. 

The  method  of  obtaining  the  electrical  connection  between 
the  ends  of  the  paper  cascade  is  shown  in  Fig.  13.   This  method 
has  the  disadvantage  that  the  operator  cannot  see  the  actual 
connection  and  must  rely  upon  the  trueness  of  the  hinges,  ma- 
chined  parts,  etc.   To  provide  sufficient  contact  pressure  a- 
long  the  joint,  C-clamps  were  used  outside  the  aluminum  bars. 

4.2  THE  CONDUCTING  MEDIUM 

Type  "L"  Teledeltos  paper  was  used  as  the  conducting  medi- 
um throughout  the  investigation.   The  specific  resistance  of 
this  paper  is  given  by  the  manufacturer  as  4000  ohms,  with  an 
accuracy  of  +.  3  to  +  8  percent  with  the  lower  value  of  resist- 
ance along  the  length  of  the  roll  and  the  higher  value  across 
the  width  of  the  roll.   The  quantity  specific  resistance,  is 
defined  and  explained  in  Appendix  D.   The  corresponding  specif- 
ic resistance  of  the  copper  bus  bars  is  clearly  negligible  in 
comparison  with  that  of  the  paper.   At  the  hinged  ends  of  the 
table,  and  near  the  electrical  joint  on  the  under  surface,  where 
the  bus  bars  are  necessarily  interrupted,  the  required  conduc- 
tivity was  secured  by  painting  the  paper  on  its  'conducting  side 

ith  a  special  silver  paint  supplied  by  the  manufacturer.   The 
specific  resistance  of  such  painted  areas  is  given  by  the  manu- 
facturer as  about  one  to  four  ohms.   While  this  value  is  consid- 
erably higher'  than  that  for  the  copper  bus  bars,  it  is  still 
negligible  in  comparison  with  that  of  the  paper.   Thus  the  con- 
ductivity of  either  the  bus  bars  or  of  the  painted  areas  is  es- 
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sentially  infinite  in  comparison  with  that  of  the  paper.   Each 
test  setup  required  a  new  piece  of  paper  making  it  a  consumable 
item  in  the  investigation. 

4.3   THE  BLADE  TEMPLATES 

The  following  profiles  were  selected  for  investigation: 

a.  The  Mutterperl  profile,  ordinates  presented  in 
Ref .  7  and  listed  in  Table  II 

b.  A  typical  compressor  profile,  the  NACA  65-( 15)10 
section,  ordinates  presented  in  Ref.  9  and  listed 
in  Table  III 

Coordinates  of  these  profiles,  scaled  to  the  chord  sizes  used 
in  this  investigation,  are  presented  in  Tables  IV  and  V. 

In  order  to  transpose  the  coordinates  to  the  paper  in  the 
proper  geometry  and  to  maintain  accuracy,  templates  were  used  as 
cutting  guides.   They  were  made  of  1/16  inch,  oil  hardening  flat 
stock  which  was  left  in  the  spheroidized  condition  instead  of 
being  hardened,  thus  avoiding  distortions. 

Tolerance  from  the  tabulated  coordinates  was  held  to  ±.005 
inches.   Each  template  was  checked  dimensionally  at  fifteen  dia- 
meters magnification  on  the  Bausch  and  Lomb  Contour  Measuring 
Projector,  Model  33-12-01,  at  the  Supersonic  Wind  Tunnel  of  the 
U.  S.  Naval  Postgraduate  School. 

To  position  the  templates  for  cutting  on  the  cascade  layout 
lines,  three  or  more  locator  holes  were  drilled  in  each.   These 
holes  were  spaced  along  a  line  parallel  to  the  chord  and  offset 
from  it.   One  of  the  holes  was  located  at  midchord.   These  pro- 
files are  shown  in  Fig.  14. 
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4.4   ELECTRICAL  CIRCUITRY  AND  INSTRUMENTS 

The  electrical  circuit  of  the  continuous  cascade  table  and  of 
the  associated  instruments  for  measuring  percentage  potential 
is  shown  in  Fig.  15.   The  power  supply  used  was  a  Heathkit  Bat- 
tery Eliminator,  Model  BE-4.   The  voltage  output  was  maintained 
at  about  11  volts  in  order  to  give  a  quicker  and  more  definite 
response  on  the  null  meter.   The  null  meter  used  was  a  General 
Electric  DC  Micro  Ammeter,  Type  DO-71,  Model  A292AA1  of  10  mi- 
croamperes full  scale  deflection.  At  a  100  percent  potential 
equal  to  about  11  volts,  C.01  percent  potential  was  the  minimum 
discernable  change  on  this  meter.   The  voltage  divider,  manu- 
factured by  Electro-Measurements  Inc.,  was  a  hexadecade  type; 
Model  RV622.   Its  linearity  as  cited  by  the  manufacturer  was 
+  0.0025  percent  with  an  input  resistance  of  10,000  ohms.   As 
mentioned  above,  0.01  percent  potential  was  the  minimum  detect- 
able change  and  therefore  only  the  four  higher  decades  of  the 
voltage  divider  were  used.   The  test  probe  or  stylus  was  sup- 
plied with  the  General  Electric  Analog  Field  Plotter  Kit.   It 
was  fitted  with  a  shallow  tapered  needle  from  a  drafting  in- 
strument . 

The  approximate  derivative  of  the  mapping  function  could 
be  determined  by  graphical  means  or  by  direct  measurement.   The 
latter  method  was  used.   To  accomplish  this  determination,  the 
potential  difference  over  a  small  distance  was  measured  by  a 
double  electrical  probe.   This  probe  was  adapted  from  a  small 
bow  compass  by  placing  a  second,  insulated  needle  in  the  lead 
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receptacle  and  soldering  both  needles  to  hookup  wires.   The  com- 
pass was  adjusted  to  the  smallest  needle  gap  possible  and  locked 
in  position.   This  gap  measured  0.04048  inches  on  the  measuring 
projector  previously  mentioned.   The  entire  compass  was  finally 
covered  with  electrical  insulating  tape. 

Fig.  16  shows  the  circuit  schematic  for  determining  the  ap- 
proximate derivative  of  the  mapping  function.   The  hookup  wires 
from  the  probe  were  lead  to  a  Leeds  and  Northrup  Portable  Pre- 
cision Potentiometer  No.  8662,  where  they  were  connected  to  the 
E.M.F.  binding  posts.   This  instrument  has  a  range  of  from  0  to 
85  millivolts.   The  power  supply  was  two  six  volt  batteries  con- 
nected in  series  to  give  approximately  12  volts.   Practically  all 
stations  on  the  profile  gave  voltage  increments  across  the  double 
probe  which  were  within  the  range  of  the  potentiometer  at  a  bat- 
tery supply  of  about  12  volts.   The  voltage  across  the  cascade 
(100  percent  potential)  was  measured  to  the  nearest  0.01  volts 
on  a  Weston  DC  Voltmeter,  Model  45,  having  0-15,  0-150  volts  full 
scale  deflection. 
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5.   TEST  PROGRAM  AND  PROCEDURE 

The  test  program  was  designed  to  show  the  effects  of  various 
geometric  parameters  upon  the  performance  of  a  cascade.   The  param- 
eters  which  define  the  geometry  of  two  dimensional,  rectilinear, 
infinite,  cascades  are  the  blade  profile,  the  stagger  angle  6, 
and  the  solidity  c/s .   The  quantities  6  ,  c  and  s  are  illustrated 
in  Fig.  17. 

5.1  DEFINITION  OF  TERMS 

The  direction  of  the  incoming  flow  may  be  specified  in  terms 
of  the  approach  angle  dC,    ,  which  is  measured  from  a  line  perpen- 
dicular to  the  cascade. 

An  alternative  parameter  which  is  sometimes  useful  is  the 
absolute  angle  of  attack  oCq  ,  which  is  defined  as  the  angle  be- 
tween the  mean  flow  direction  and  the  direction  of  zero  lift, 
Axis  I .   The  mean  flow  angle  oCm   »  depends  both  upon  the  inflow 
and  outflow  angles  c^  and  cCz    .   The  angle  of  zero  lift  dC{jQ   , 
of  the  cascade  is  not  initially  known,  since  it  depends  on  the 
entire  cascade  geometry  in  some  complex  manner,  but  it  may  be 
found  from  the  results  of  the  test  measurements  in  any  given 
case.   Once  cCu       has  been  determined,  a  convenient  parameter  for 
further  analysis  of  the  cascade  is  the  relative  angle  of  zero 
lift  p  ,  which  is  defined  as  the  angle  between  the  chord  and 

Axis  1,  the  direction  of  zero  lift  of  the  profile  in  cascade. 

These  various  relations  are  shown  in  Fig.  17. 

The  performance  of  a  cascade  may  be  expressed  in  terms  of 
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any  convenient  parameter  which  relates  the  inflow  and  outflow 
directions.   Among  these  are  the  turning  angle  ©     and  the  lift 
coefficients  /C^  ,  >o   or  /Zo      •   The  turning  angle  S    ,    repre- 
sents  the  net  change  in  direction  of  flow  from  inlet  to  exit. 
The  lift  coefficients  C^    ,  *,£  and  /C^        represent  the  ra- 

tio of  the  mean  effective  pressure  difference  across  the  blade 
to  some  reference  dynamic  pressure.   They  differ  only  in  the  ref- 
erence condition  used:  /Co     uses  the  mean  velocity  head,  while 
/Co      and  &g        employ  the  inlet  and  outlet  velocity  heads  re- 
spectively.  All  three  of  these  coefficients  are  closely  relat- 
ed to  the  turning  angle,  and  the  choice  of  parameters  is  large- 
ly a  matter  of  convenience,  although  for  a  real  fluid  the  quan- 
tity £,p       appears  to  be  an  approximate  criterion  of  blade 
stall. 

It  can  be  shown  theoretically  that  for  an  isolated  airfoil 

/C%   =  ZTt  sanofa  5.1 

Similarily,  for  a  cascade  of  airfoils,  it  can  be  shown  that 

/Ci  =  ZrrKstnoCa 

where  oCq   =  o^n  —  o^  and  K.  is  the  cascade  coefficient  which  is  a 
constant  determined  by  cascade  geometry. 

There  are  two  basic  parameters,  then,  which  completely  de- 
termine the  performance  of  a  given  cascade  in  ideal  fluid  flow. 
These  are  K  and   P)    .   They  have  been  calculated  for  the  cas- 
cades tested  and  are  reported. 
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5.2   THE  TEST  PROGRAM 

The  method  used  in  this  investigation  requires  that  a  sep- 
arate set  of  measurements  be  made  for  each  cascade  configuration 
tested,  whereupon  the  effect  of  variations  in  inflow  direction 
can  be  determined  analytically.   This  greatly  reduces  the  amount 
of  testing  that  is  required;  nevertheless,  the  testing  procedure 
was  necessarily  painstaking,  and  due  to  the  developmental  nature 
of  this  investigation  only  a  limited  number  of  tests  could  be 
undertaken  and  analyzed  in  the  time  available  for  this  study. 

The  various  combinations  of  blade-  profile,  solidity  and 
stagger  angle  which  were  investigated  are  summarized  in  Table  VI 
and  discussed  below. 

The  airfoil  shape  presented  by  Mutterperl  in  Ref .  7  was  se- 
lected for  the  first  series  of  tests.   Published  data  on  this 
profile  had  been  determined  analytically.   Tests  on  the  profile 
in  this  investigation  included  solidities  of  1.0  and  1.5,  both 
at  a  stagger  angle  of  45°.   The  latter  case  was  not  covered  by 
Mutterperl  in  his  report. 

The  second  series  of  tests  was  designed  around  a  typical 
compressor  profile  at  combinations  of  stagger  angle  and  solidity 
characteristic  of  practical  compressor  geometry  (based  on  the 
authors'  experience  with  the  Gas  Turbine  Division  of  General 
Electric  during  the  1957  Summer  Field  Trip  from  the  U.  S.  Naval 
Postgraduate  School). 

Two  groups  of  tests  were  conducted  using  an  NACA  65-( 15)10 
profile  at  stagger  angles  of  20°,  25°,  30°,  and  35°.  One  group 
used  the  six  inch  chord  at  a  solidity  (T   -    1.0  while  the  other 
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used  the  twelve  inch  chord  at  (T  =  1.5.   At  least  one  check 
point  for  each  test  was  available  from  the  experimental  data 
presented  in  Ref .  9. 

One  of  the  objectives  of  this  investigation  was  to  demon- 
strate an  extension  of  the  method  namely,  the  determination  of 
the  velocity  distribution  around  a  blade  in  cascade.   Data  for 
this  demonstration  was  gathered  during  the  test  of  NACA  65- (15) 
10  (12"  chord)  at  (T  =  1.5  and  Q    =  35°.   Selection  of  this  test 
provided  the  opportunity  to  compare  results  with  those  publish- 
ed in  Ref.  9. 

5 . 3  PROCEDURE 

The  procedure  is  perhaps  best  illustrated  by  following  a 
typical  cascade  setup  through  to  data  collection. 

The  first  step  of  the  procedure  was  preparation  of  the  con- 
ducting paper  for  the  table.   The  paper  was  supplied  in  rolls, 
34  inches  wide  and  about  20  feet  long.   The  working  surface  had 
been  designed  with  the  electrodes  positioned  to  accommodate 
this  width  so  it  was  necessary  only  to  trim  the  paper  to  length. 
The  length  required  varied  with  solidity  for  a  given  blade  chord. 
For  (T  =  1.0  and  a  six  inch  chord,  88  inches  of  paper  were  used. 
For  the  tests  with   (T  =  1.5  and  the  twelve  inch  chord,  92  in- 
ches were  required. 

The  paper  was  laid  out  flat  on  a  smooth  table  and  fastened 
without  stretching  by  masking  tape  at  intervals  along  the  edges. 
The  black  (non-conducting)  surface  was  exposed  as  the  drawing 
surface  for  layout  lines. 
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Fig.  18  shows  a  sketch  of  the  layout  lines  on  a  typical 
piece  of  paper.   Fig.  19  shows  the  template  outline  on  one  set 
of  these  lines.   Included  in  this  latter  sketch  are  the  refer- 
ence lines  near  the  leading  and  trailing  edges  which  were  used 
to  locate  data  points. 

The  cascade  center  line  was  drawn  first.   The  best  results 
in  keeping  this  line  straight  occured  when  using  a  chalkline  of 
sewing  thread  and  blackboard  chalk.   The  line  so  obtained  was 
then  penciled  in  using  a  six  foot  straight  edge.   Exact  center- 
ing of  this  line  on  the  paper  was  not  necessary  as  the  copper 
electrodes,  being  32  inches  apart,  allowed  some  deviation. 

The  cascade  center  line  was  used  as  the  reference  for  lay- 
out of  the  selected  spacing  and  stagger  angles.  Care  was  taken 
to  locate  the  profiles  so  that  the  electrical  joint  would  occur 
midway  between  adjacent  blades  along  this  reference  line.  In- 
stances of  large  stagger  angle,  where  the  seam  passed  through 
corresponding  sections  of  the  cutouts  on  each  end  of  the  paper, 
could  be  handled  by  the  joining  device. 

A  drafting  machine  was  utilized  to  construct  the  angular 
stagger  and  midchord  lines  shown  in  Figs.  18  and  19. 

In  order  to  maintain  accuracy  in  locating  data  points,  ref- 
erence lines  were  drawn  on  the  three  center  profiles  of  the  pa- 
per cascade.   These  lines  were  constructed  perpendicular  to  the 
chord  and  extended  to  intersect  the  contour  of  the  profile  cut- 
outs as  shown  in  Fig.  19.   The  interval  along  the  chord  of  these 
lines  was  1/40  inch  for  the  six  inch  chord  and  1/20  inch  for  the 
twelve  inch  chord.   A  jewelers'  eyepiece,  and  6H  pencil,  and  the 
scale  of  the  drafting  machine  arm  were  used  to  construct  these 
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lines.   Only  part  of  the  profile  was  lined  in  this  manner.   Fif- 
teen lines  were  drawn  at  the  trailing  edge  in  all  tests.   The  num- 
ber and  position  of  the  lines  near  the  leading  edge  varied  with 
the  location  of  the  front  stagnation  point. 

The  paper  was  then  removed  to  a  Masonite  surface  for  cutting 
out  the  profiles.   In  this  operation  the  template  was  located  by 
means  of  the  holes  provided.   Cutting  proceeded  from  the  tail 
and  from  the  nose  toward  the  midchord  where  the  cuts  met.   The 
knife  used  for  cutting  was  an  "Xacto"  using  the  blade  having  the 
straight  cutting  edge.   The  blade  was  honed  to  a  shape  illus- 
trated in  Fig.  20  and  stropped  on  a  Masonite  surface  to  remove 
all  burs.   In  cutting,  it  was  imperative  to  hold  the  knife  so 
the  cutting  edge  was  at  an  angle  of  about  45  to  the  paper  and 
to  guide  it  in  toward  the  template.   All  other  methods  of  hold- 
ing and  cutting  either  did  not  adhere  to  the  profile  shape  or 
produced  jagged  edges. 

The  ends  of  the  paper  were  then  folded  to  fit  over  the 
corner  of  the  aluminum  bars  at  the  joint  of  the  table.   The  fold 
line  was  marked  on  the  non-conducting  (black)  surface  so  the  pa- 
per was  first  folded  to  keep  this  line  exposed  as  a  guide.   After 
this  fold  had  been  sharply  creased,  it  was  reversed  to  expose  the 
conducting  surface. 

Painting  the  bus  bar  extensions  was  the  next  step  in  the 
preparation  of  the  paper  cascade.   To  provide  guidelines  for  this 
painting,  the  paper  was  marked  while  in  its  proper  position  on 
the  working  surface.   A  center  line  had  been  scribed  on  the  table 
during  its  construction.   Bringing  the  cascade  center  line  into 
coincidence  with  the  scribed  line  properly  positioned  the  paper 
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in  one  direction.   Positioning  in  the  other  direction  was  ac- 
complished by  fitting  the  folds  at  each  end  of  the  paper  to  the 
right  angle  corner  on  their  respective  parts  of  the  joining  de- 
vice.  This  is  shown  in  cross  section  in  Fig.  13,  with  the  paper 
joined.   The  end  sections  of  the  paper  were  held  by  suitable 
weights  on  each  end  section  of  the  table  while  the  middle  of  the 
paper  was  moved  along  the  scribe  line  of  the  main  table  s'ection 
until  the  bights  of  the  paper  over  the  hinges  were  equalized. 
Fig.  12  shows  the  paper  clamped  in  this  position.   The  guides 
for  the  bus  bar  extensions  were  then  marked  on  the  paper's  con- 
ducting surface  by  pressing  it  down  over  the  sharp  inner  edge  of 
the  ends  of  the. copper  bus  bars. 

The  paper  was  then  removed  to  the  drawing  surface  and  layed 
out  flat  with  its  conducting  surface  up.   At  each  end  of  the  pa- 
per the  painting  guides  were  pencil  lines  which  had  been  extend- 
ed along  the  bus  bar  impressions  to  the  ends  of  the  paper. 
Painting  guides  for  the  bights  of  the  paper  were  obtained  by  pen- 
cil lines  connecting  the  bus  bar  impressions  at  the  ends  of  each 
bight.   Scotch  Brand  Electric  Tape  was  used  as  a  masking  tape. 
It  was  placed  along  the  above  lines  on  the  side  away  from  the 
paper  edge  and  extended  beyond  the  bus  bar  marks.   The  edges  of 
the  paper  were  then  painted  with  the  silver  conducting  paint  a- 
long  each  of  the  masked  areas.   When  the  paint  had  set  but  not 
dried,  the  masking  tape  was  removed.   Unless  this  tape  had  been 
dusted  lightly  with  blackboard  chalk  or  some  other  dust,  the 
graphite  conducting  medium  was  liable  to  pull  away  from  the  pa- 
per during  tape  removal. 

At  this  point  various  reinforcing  sections  of  tape  were  ap- 
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plied  to  the  paper.   This  was  especially  needed  when  using  the 
twelve  inch  profile  cutouts. 

After  the  paint  had  dried,  the  paper  was  aligned  on  the  ex- 
tended table,  non-conducting  surface  up  as  when  marking,  and 
held.   The  electrical  tape  was  used  in  small  pieces  to  hold  the 
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paper  in  position  by  placing  it  along  the  paper  edges  behind  the 
contact  line  of  the  electrodes  and  over  blade  cutouts  not  used 
for  data.   The  clamps  were  then  applied  as  shown  in  Figs.  10  and 
12. 

The  ends  of  the  paper  sticking  through  the  slots  of  the  a- 
luminum  bars  were  cut  off  to  leave  about  3/8  inch  of  paper  ex- 
tending beyond  the  fold.   A  3/8  inch  wide  strip  of  electrical 
tape  applied  to  the  cut  edge  of  these  ends  was  used  to  hold  the 
fold  over  the  corner  of  its  respective  bar.   The  relief  in  the 
movable  bar  gave  room  for  the  tape  and  at  the  joint  itself  only 
the  two  layers  of  pa'per,  conducting  sides  together,  remained. 

The  table  was  turned  over  and  the  ends  folded  toward  each 
other.   The  folds  of  the  paper  were  joined  simultaneously  and  a 
support  under  the  two  ends  prevented  their  moving  through  their 
closest  point  of  approach.   The  ends  of  the  table  were  secured 
against  the  middle  by  bolts.   Insulating  paper  was  then  slipped 
between  the  conducting  paper  and  aluminum  bar  and  pushed  to  the 
joint  to  prevent  the  curled  cutout  edges  from  contacting  the  a- 
luminum  bars.   Oclamps  were  applied  to  press  the  bars  together 
and  make  a  tight  contact  between  the  folds  on  the  ends  of  the 
paper. 

The  table  was  turned  over  and  the  circuit  connected  as  in 
Fig.  15. 
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The  data  points  were  numbered  from  leading  edge  to  trailing 
edge.   Stations  every  1/8  inch  along  the  chord  on  the  six  inch 
profile  were  numbered  consecutively.   Data  points  between  these 
lines  were  lettered  a,  b,  c,  and  d.   The  data  consisted  of  the 
value  of  O   at  the  leading  edge,  the  points  near  Oi  »  and  at 
the  fifteen  upper  and  lower  points  near  the  trailing  edge.   To 
obtain  ^       ,  a  line  perpendicular  to  the  electrodes  and  the 
cascade  center  line  was  drawn  from  the  electrode  ahead  of  the 
cascade.   It  was  marked  at  intervals  of  one  inch  for  five  or  six 
inches  from  the  bus  bar.   The  value  of  O   was  recorded  for  each 
of  these  points.   All  values  of  (S   were  read  and  recorded  to 
0.C1  percent  potential. 

When  obtaining  the  derivative  approximation  ~-  it  was 
necessary  to  use  the  circuit  of  Fig.  16.   Percent  potential 
could  be  obtained  simultaneously  by  connecting  the  voltage  di- 
vider with  null  meter  and  probe  across  the  battery  supply. 

The  data  in  this  test  consisted  of  \J     at  various  points  a- 
round  the  profile  edge,  corresponding  values  of  AC.  in  millivolts 
obtained  from  the  precision  potentiometer  and  total  potential  E* 
in  volts.   The  double  probe  was  positioned  so  that  its  interval 
between  probes  was  centered  about  the  point  where  £>   was  obtain- 
ed and  with  its  needles  just  making  contact  at  the  edge  of  the 
profile.   The  value  of  Au.  was  read  and  recorded  to  the  nearest 
.01  millivolt. 
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6.   DEVELOPMENT  OF  RESULTS 

A  sample  of  the  data  taken  during  these  tests  for  the  deter- 
mination of  turning  angles  is  presented  in  Table  VII.   A  sample 
of  the  data  taken  for  the  determination  of  velocity  is  presented 
in  Table  VIII.   Table  IX  contains  the  input  quantities  used  for 
the  determination  of  turning  angles  herein  reported. 

Results  of  tests  on  the  cascade  analyzed  in  Ref .  7  are  pre- 
sented in  Fig.  21  as  a  curve  of  coefficient  of  lift  versus  mean 
velocity  direction. 

Results  of  tests  made  on  the  NACA  65-( 15)10  airfoil  are  pre- 
sented in  Figs.  22  through  31.   In  Figs.  22  and  25  the  values  of 
turning  angle  for  various  blade  stagger  angles  are  plotted  versus 
cascade  entrance  angle  with  lines  of  constant  coefficient  of  lift 
/C/   crossplotted.   Fig.  22  is  for  a  solidity  of  1.0  and  Fig.  25 
is  for  a  solidity  of  1.5.   Figs.  23  and  26  contain,  for  a  solid- 
ity of  1.0  and  1.5  respectively,  curves  of  coefficient  of  lift  /C^ 
for  various  blade  stagger  angles,  as  a  function  of  cascade  en- 
trance angle.   Figs.  24  and  27  show,  for  a  solidity  of  1.0  and  1.5 
respectively,  curves  of  coefficient  of  lift  C%       for  various  blade 
stagger  angles,  as  a  function  of  cascade  entrance  angle.   Fig.  28 
is  a  plot  of  cascade  coefficient  K  versus  blade  stagger  angle  £> 
Fig.  29  is  a  plot  of  relative  zero  lift  angle  A    versus  blade 
stagger  angle  P>  . 

Figs.  30  and  31  contain  the  results  of  the  test  of  this  meth- 
od to  predict  velocity  distribution  about  the  NACA  65-( 15)10  air- 
foil in  cascade  at  stagger  angles  of  35°  and  solidity  of  1.5.   The 
results  are  plotted  as  [-$-)    versus  percentage  distance  along  the 
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chord . 

The  recorded  data  necessary  to  determine  turning  angles 
and  coefficients  of  lift  for  a  given  cascade  (blade,  stagger, 
and  solidity)  consisted  in  part  of  readings  of  the  percentage 
voltage  at  several  stations  of  known  location  in  the  vicinity 
of  O,  ,  Uz  ,  and  Oyr  .   See  Table  VII.   Data  were  recorded  for 
at  least  three  blades  in  each  cascade  test,  but  final  values  of 
this  report  are  based  on  the  median  blade  for  each  case.   These 
percentage  readings  near  O,  and  ^j2     could  be  plotted  versus 
distance  along  the  chord  to  permit  utilization  of  several  volt- 
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age  readings  to  insure  greater  accuracy  in  the  determination  of 
the  values  of  £//  and  o2  .   See  Fig.  32. 

Finding  Q>yi   vvas  a  more  difficult  problem.   Fig.  3  indicates 
that  the  curve  of  O  versus  percentage  of  the  chord  becomes  tan- 
gent to  the  line  representing  100%  chord  at  the  trailing  edge. 
Determination  of  O^.  ,  the  point  of  tangency,  by  experimental 
measurement  meant  reliance  upon  a  single  potential  reading  be- 
cause of  the  physical  nature  of  the  very  sharp  trailing  edge. 

The  method  used  in  this  investigation  to  determine  O7*  more 
accurately  involved  obtaining  an  analytical  curve  of  C   near  the 
trailing  edge  of  the  airfoil.   The  form  of  the  curve  was  that  of 
potential  flow  around  a  straight  sided  wedge.   The  curve  was  ob- 
tained by  fitting  this,  form  to  several  experimental  points  in  a 
finite  but  small  distance  near  the  trailing  edge.   Curve  fitting 
was  accomplished  by  means  of  the  method  of  least  squares.   The 
point  of  tangency  C,jr    was  then  determined  analytically.   Appen- 
dix E  contains  the  details  of  this  curve  fitting  procedure. 

The  value  of  °^      was  obtained  by  plotting.   Near  the  bus 
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bars  readings  of  o  were  plotted  versus  perpendicular  distance 

from  the  bus  bars.   The  slope  of  this  straight  line  was  ^v 

See  Fig.  32. 

The  determination  of  the  velocity  at  a  point  on  the  blade 

required  the  measurement  of  AL.  at  that  point.   A  sample  of 

recorded  readings  of  £>  fc.  is  contained  in  Table  VIII.   These 

readings  were  converted  to 


L*  d 


by  means  of  the  relation 
where 


d  is  the  separation  of  the  double  probe.   Values  of 


A& 
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are 


also  shown  in  Table  VIII.   Values  of 
by  means  of  the  formula 


V. 


were  then  obtained 


V?  -  /X/  COS  oC, 
V,    A?  VQ 


1   ~mW^to nofm  +  "rfri$r&- 


6.2 


The  NCR  102A  electronic,  digital  computer  of  the  U.  S.  Na- 
val Postgraduate  School  was  used  to  compute  o<t»  , oC\  ,  oCi  »  oCi     » 
and  -Q-    .   Machine  assistance  proved  indispensible  for  the  ex- 
tensive  calculations  required  by  the  many  cascade  combinations 
tested  during  this  investigation. 

One  program  on  the  computer  was  used  to  determine  Of   •   In- 
put data  consisted  of  twenty  potential  readings,  ten  from  the 
upper  surface  on  the  airfoil  near  the  trailing  edge  and  ten  from 
the  lower  surface  near  the  trailing  edge,  together  with  their 
respective  distances  from  the  trailing  edge  and  the  value  of  p. 
See  Appendix  E.   The  computer  solved  for  Gj,  . 

A  second  computer  program  was  used  to  determine  oCt    ,  qCz    , 
and  d(i     -      Input  data  consisted  of  u,  ,  (jA  ,  (J,    ,  {J      ,   5  ,  and 
f-^p-j  •   in  addition,  inputs  of  O   ar>d  ^%-  were  required  for 
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each  of  the  points  at  which  velocity  ratios 


were  desired. 
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7.      DISCUSSION 

7.1  GENERAL 

The  discussion  of  this  method  for  predicting  the  flow  of  an 
ideal  fluid  across  any  arbitrary  cascade  first  examines  the  re- 
sults of  the  tests  made  on  various  cascades.   The  analysis  of 
the  results  then  provides  the  basis  for  a  critique  of  the  method 
and  for  suggested  improvements  in  equipment  and  in  technique. 
Possible  applications  of  this  method  to  other  aspects  of  engi- 
neering are  also  suggested. 

7 . 2  RESULTS 

The  geometric  parameters  which  determine  a  cascade  are  the 
blade  profile,  the  stagger  angle,  and  the  solidity.   For  a  given 
cascade  the  quantities  K  and  yu     are  the  fundamental  parameters 
which  were  determined  by  this  investigation  for  indicating  the 
performance  of  a  cascade  over  a  wide  range  of  mean  velocity  di- 
rections.  Turning  angles  were  also  calculated  for  each  cascade. 
From  this  information  other  parameters'  were  calculated  to 

4 

permit  cascade  analysis  from  different  viewpoints  and  to  permit 
comparison  with  published  reports. 

It  first  seemed  appropriate  to  compare  the  results  of  the 
method  herein  discussed  with  a  published  analytical  solution. 
Ref.  7  contains  an  analytical  solution  of  potential  flow  across 
a  given  cascade  of  (T  =  1  and  P  =  4b.   This  report  was  used  for 
comparative  purposes  because  it  contained  coordinates  for  the 
profile  used  and  because  results  were  reported  over  a  range  of 
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directions  of  the  mean  velocity.   The  blade  profile  used  is  sim- 
ilar in  camber  and  thickness  to  those  presently  found  in  compres- 
sors.  The  blade,  however,  has  a  trailing  edge  of  zero  thickness. 
The  pointed  trailing  edge  is  difficult  to  duplicate  exactly  in 
the  paper  cutouts  and  would  not  be  used  in  actual  compressors. 
Reference  to  Fig.  21  indicates  that  the  method  of  this  paper  ap- 
parently predicts  coefficients  of  lift  which  are  about  14%  great- 
er than  those  contained  in  Ref.  7.   The  difference  represents  a- 
bout  5°  in  turning  angle.   This  difference  was  greater  than  had 
been  expected,  but  it  is  subsequently  shown  that  this  disparity 
can  be  reduced.   Subsequent  tests  indicated  that  an  increased 
error  in  turning  angle  resulted  from  tests  on  cascades  with  large 
stagger  angles.   The  larger  stagger  angle  cascades  are  more  sus- 
ceptible to  errors  caused  by  the  variation  of  specific  resistance 
of  the  paper  between  width  and  length  directions.   Possible  solu- 
tions to  this  problem  are  discussed  later. 

Fig.  21  does  indicate  that  the  slopes  of  the  curves  of 
C£     vs .  oC     are  almost  the  same.   Calculations  reveal  that  the 
cascade  of  Ref.  7  has  a  coefficient  K  of  .74  while  the  cascade 
coefficient  determined  by  this  method  is  .72.   Thus  if  determi- 
nation of  the  angle  of  zero  lift  were  made  separately,  the  turn- 
ing angle  predicted  by  this  method  would  be  within  a  fraction  of 
a  degree  of  the  turning  angles  reported  in  Ref.  7. 

Fig.  21  indicates  that  increase  of  the  solidity  of  this 
cascade  to  1.5  changes  the  cascade  coefficient  to  .49.   The  1.5 
solidity  cascade  is  not  analyzed  in  Ref.  7. 

Ref.  9  is  a  report  of  systematic  two  dimensional  cascade 
tests  of  NACA  65  series  compressor  blades.   Cascades  considered 
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to  be  most  typical  of  present  compressor  design  were  selected 
from  the  report  for  investigation  by  this  method. 

The  NACA  test  data  was  taken  by  varying  cascade  stagger  an- 
gle. This  was  required  because  of  the  difficulty  of  varying  the 
direction  of  the  approaching  velocity  in  the  two  dimensional 
wind  tunnel.  The  method  of  this  report  permits  variation  of  the 
direction  of  approaching  velocity  for  constant  cascade  geometry. 
Comparison,  then,  was  possible  only  at  certain  points,  indicated 
in  Figs.  22,  23,  25  and  26. 

It  must  be  borne  in  mind  that  Ref .  9  is  a  report  of  tests 
performed  using  an  actual  fluid.  The  magnitude  of  the  effect  of 
viscosity  is  uncertain,  but  one  could  not  expect  results  based 
on  viscous  fluid  flow  to  be  identical  to  those  for  ideal  fluid 
flow.  Comparison  serves  to  indicate,  nevertheless,  that  pre- 
dictions based  upon  ideal  fluid  flow  are  sufficiently  close  to 
measured  values  for  air  to  be  of  value  in  actual  cascade  design. 

Figs.  22  and  25  give  the  results  for  the  cascade  of  the 
NACA  65-(15)10  airfoils  with  (T  =  1  and  1.5  respectively.   The 
figures  show  that  the  relationship  between  turning  angle  and  cas- 
cade entrance  angle  is  almost  a  linear  one,  indicating  that  the 
exit  angle  remains  nearly  constant.   For  a  given  cascade  en- 
trance angle,  turning  angle  is  of  course  increased  by  decreasing 
the  stagger  angle.   The  effect  of  increasing  solidity  on  turning 
angle  (and  therefore  on  the  quantity  (T/Cp  )  is  small.   Hence  the 
cross  plots  indicate  an  appreciable  reduction  in    /Co  with 

increase  in  solidity.   These  are  all  general  results  which  one 
would  expect. 

The  cascade  designer  wonders  what  represents  optimum  condi- 
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tions  and  the  extent  of  variation  permissible  for  satisfactory 
performance.   The  cross  plotted  values  of  Co        give  an  indica- 
tion of  the  coefficient  of  lift  associated  with  a  certain  turn- 
ing angle.   This  information  together  with  the  plot  of  ideal  an- 
gle of  attack  helps  to  denote  an  area  of  satisfactory  operation 
based  upon  upstream  conditions.   However,  for  the  decelerating 
flow  experienced  in  a  compressor,  information  showing  the  extent 
of  the  range  of  operation  for  an  optimum  drag-lift  ratio  is  de- 
sired.  For  this  reason  Figs.  24  and  27  have  been  presented. 
Ref.  11  says  tests  show  that  the  drag-lift  ratios  of  cascades 
varying  from  turbine  to  compressor  types  show  a  minimum  at  an 
approximately  constant  value  of  /C^   of  about  .8.   For  blades 
of  the  types  generally  used  in  cascades  the  drag  coefficient  is 
low  and  almost  constant  over  a  range  which  extends  perhaps  five 
degrees  above  and  fifteen  degrees  below  the  minimum  drag-lift 
ratio.   These  general  criteria  for  determining  range  of  opera- 
tion are  indicated  on  Figs.  22  and  25. 

The  effect  of  changing  cascade  stagger  is  apparent  from 
Figs.  22  and  25.   One  would  expect  the  spacing  between  lines  of 
constant  &     to  be  more  uniform  than  this  test  indicated.   This 
variation  is  considered  to  be  a  manifestation  of  the  experimen- 
tal error  of  the  method.   For  the  whole  series  of  tests  the  error 
is  believed  to  be  +  1°.   A  test  program  which  included  a  wider 
range  of  stagger  angles  or  which  provided  for  smaller  test  in- 
crements in  this  same  range  would  serve  to  assist  in  reducing 
this  random  error  by  permitting  the  application  of  statistical 
data  smoothing  techniques. 

There  is  again  an  indication,  particularly  for  the  cascade 
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test  at   (T  =  1.5  and   u     -   35°  that  the  high  stagger  angle  cas- 
cade is  subject  to  an  additional  error.   The  turning  angles  pre- 
dicted for  this  case  seem  to  be  greater  than  one  would  expect, 
considering  the  other  three  cascades.   As  in  the  45°  stagger  an- 
gle cascade  first  discussed,  some  of  this  difference  is  probably 
attributable  to  the  paper. 

This  particular  cascade  was  also  analyzed  to  see  what  value 
of  O-tt  would  give  results  in  agreement  with  the  NACA  report. 
The  point  on  the  blade  cutout  which  had  the  value  of  U      so  de- 
termined was  located  on  the  upper  side  of  the  blade  profile  very 
near  the  trailing  edge.   This  suggests  that  the  real  fluid  in 
the  NACA  tests,  under  the  influence  of  viscosity  and  adverse 
pressure  gradient,  separated  from  the  airfoil  at  a  point  closely 
approximating  this  same  location  on  the  blade  cutout. 

For  the   (T  =  1  cascade  the  predicted  turning  angles  com- 
pare remarkably  well  with  the  NACA  test  points  for  a  cascade  en- 
trance angle  c(,   -   45.   The  comparison  is  not  as  close  for  the 
two  NACA  points  plotted  at  a  cascade  approach  angle  cC,   =  30°. 
For  y     =  20°  a  line  drawn  between  the  two  NACA  test  points  at 
ocj  =  30  and  oCi  =  45°  would  not  have  the  same  slope  as  the 
corresponding  line  produced  by  this  method.   The  same  thing  is 
true  to  a  lesser  extent  for  p  =  25°.   These  differences  in 
slope  are  considered  to  be  a  consequence  of  viscosity. 

For  the   (T  =  1.5  cascade  there  is  a  somewhat  larger  differ- 
ence between  predicted  turning  angles  and  reported  NACA  test  re- 
sults, but  the  comparison  is  still  close.   Twelve  inch  chord 
blade  templates  were  used  for  three  of  the  four  cascades  of 
vT  =  1.5.   Only  the   p  =  25°  cascade  was  cut  out  using  a  six 
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inch  template.   The  somewhat  increased  error  due  to  bus  bar 
location  for  the  larger  template  is  not  discernable;  so  use  of 
the  larger  blade  is  probably  warranted  for  high  solidity  cas- 
cades when  the  reduction  in  cut  out  time  and  in  labor  are  con- 
sidered.  A  series  of  analyses  for  various  values  of  p  and  (T 
such  as  is  contained  in  Appendix  F  for   (T  =  1  and   p=  45° 
would  permit  approximate  determination  of  the  largest  sized 
blade  permissable  for  a  given  acceptable  percentage  error. 

Fig.  28  presents  a  concise  summary  of  the  effects  of  solid- 
ity and  stagger  angle  upon  cascade  performance.   The  variation 
of  cascade  coefficient  K  with  changes  in  stagger  angle  are 
shown  for  each  solidity.   For  purposes  of  approximate  comparison 
values  of  K  for  a  flat  plate  cascade  having  the  same  nominal 
value  of  solidity  and  stagger  as-  the  actual  cascade  are  also 
plotted  and  are  seen  to  be  quite  similar. 

Another  concise  representation  of  cascade  characteristics 
is  contained  in  Fig.  29  which  shows  the  effect  of  changing  stag- 
ger angle  upon  the  relative  angle  of  zero  lift  oC^     .   This  pre- 
sentation amplifies  and  clearly  shows  the  variation  in  the  an- 
gle of  zero  lift  with  changes  in  stagger  angle.   It  suggests  that 
a  similar  plot  for  a  much  larger  number  of  tests  would  permit  the 
determination  of  a  smooth  curve  of  p   vs . 6 .   If  the  magnitude 
of  random  error  were  thus  partially  removed  from  the  determina- 
tion of  o^    ,  a  more  uniform  spacing  of  the  constant  ft  lines 
in  Figs.  22  and  25  would  result  and  the  error  indicated  in  Fig. 
21  would  be  considerably  reduced.   Unfortunately  time  did  not 
permit  conducting  the  additional  number  of  tests  which  would  be 
required  to  determine  this  curve. 
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The  curve  of  Fig.  31,  rrf~~)     vs.  percent  chord  compares  favor- 
ably with  the  NACA  curve  for  a  similar  cascade  entrance  angle. 
However,  Fig.  31  is  a  smooth  curve  determined  from  calculations 
based  upon  raw  data  by  means  of  a  method  of  least  squares  de- 
veloped by  Professor  Gawain.   Before  extensive  tests  of  veloc- 
ity could  be  undertaken,  the  double  probe  and  perhaps  the  tech- 
nique of  'obtaining  &L-   readings  would  require  improvement.   Fig. 
31  does  indicate,  however  that  the  velocity  distribution  can  be 
determined  by  this  method. 

7 . 3  ACCURACY 

The  accuracy  of  the  results  obtained  by  this  method  of 
cascade  analysis  is  a  function  of  the  approximations  made  and 
of  the  accuracy  of  the  measured  input  data.   The  results  are 
more  sensitive  to  errors  in  some  input  quantities  than  to  errors 
in  others,  and  certain  measured  quantities  can  be  determined 
more  accurately  than  others. 

Table  X  gives  an  indication  of  the  sensitivity  of  results 
to  variation  of  input  data.   It  lists  the  change  in  calculated 
turning  angle  caused  by  individually  increasing  each  of  the  in- 
put quantities  o,  ,  Oi  »  Cf  >  an<^   4~~  by  one  percent  of  its 
measured  value.   This  information  was  compiled  for  a  cascade  of 
the  airfoil  contained  in  Ref .  7  with   6  =  45°  and   (T  =  1,  but 
it  is  considered  to  be  representative  of  the  other  blades  and 
cascades  tested.   The  table  indicates  that  the  results  are  par- 
ticularly sensitive  to  changes  of  G2  and  Of  .   It  indicates 
the  desirability  of  obtaining  the  input  measured  values  to  an 
accuracy  of  one  part  in  ten  thousand. 
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Many  operations  were  required  for  each  cascade  layout,  and 
every  step  had  some  effect  on  the  value  of  the  final  measured 
quantities.   The  largest  of  the  construction  errors  was  due  to 
the  conducting  medium.   The  value  of  the  specific  resistance  of 
this  material  was  given  by  the  manufacturer  as  400C  ohms.   This 
value  was  specified  accurate  from  +.3  to  +.8  percent  with  the 
lower  value  of  resistance  along  the  length  of  the  roll  and  the  ■ 
higher  value  across  the  width  of  the  roll.   This  variation  is 
explained  by  the  manner  of  manufacturing.   These  errors  were 
necessarily  accepted  as  the  particular  conducting  medium  used 
was  basic  to  this  investigation. 

Wood  was  used  for  structural  members  of  the  working  sur- 
face.  The  position  of  the  bus  bars  was  held  to  the  tolerance 
mentioned  in  the  equipment  section,  but  the  trueness  of  the 
joint  depended  upon  the  ability  of  the  wood  to  maintain  its  di- 
mensions under  continued  use.   It  was  found  that  after  three  or 
four  tests,  the  table  required  cleaning  because  the  graphite 
deposit  on  the  table  caused  a  change  in  the  conductivity  of  sub- 
sequent papers  when  not  removed. 

In  the  layout  of  the  cascade  on  the  paper  several  sources 
of  probable  error  were  encountered,  such  as  dimensioning,  con- 
structing layout  and  reference  lines  and  setting  stagger  angles. 
These  were  subject  to  the  same  limits  of  error  encountered  in 
engineering  drawings   A  reading  glass  and  a  jewelers'  eyepiece 
were  used  to  provide  magnification  for  greater  accuracy  in  lay- 

0 

out.   Dimensions  could  be  held  to  a  tolerance  of  .01  inch. 
Lines  were  maintained  to  widths  of  from  four  to  six  thousandths 
of  an  inch  and  angles  were  maintained  correct  to  10  minutes  of 
arc . 
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The  templates  introduced  a  systematic  error  due  to  the 
tolerance  of  manufacture  and  a  random  error  in  the  process  of 
positioning.   The  manufacturing  tolerance  was  +  0.005  inches. 
The  profile  cutouts  were  maintained -to  their  respective  coordi- 
nates with  a  tolerance  of  about  .010  inches. 

Reading  and  determination  errors  of  the  instruments  were 
encountered  only  with  the  null  meter  and  with  the  probe.   Since 
a  zero  reading  or  null  was  required  on  the  meter  prior  to  read- 
ing the  voltage  divider,  detecting  a  null  was  reduced  to  detect- 
ing motion  of  the  meter  needle.   The  sensitivity  of  the  instru- 
ments permitted  reading  potential  on  the  voltage  divider  to  one 
part  in  ten  thousand.   Readout  was  directly  in  percent  poten- 
tial; so  constancy  of  imposed  voltage  was  not  required.   Errors 
in  probe  location  depended  upon  all  previous  dimensioning  and 
construction  of  reference  lines.   With  optical  aids  it  was  pos- 
sible to  place  the  probe  within  the  actual  width  of  the  pencil 
lines.   Distance  from  the  cutout  edges  of  the  profile  could  not 
be  determined  to  comparable  accuracy.   However,  the  constant 
potential  lines  are  normal  to  the  cutout  near  the  profile  edge; 
so  this  distance  is  not  as  critical. 

The  velocity  distribution  around  the  profile  was  determined 
to  demonstrate  in  preliminary  fashion  one  of  the  capabilities  of 
the  method,   this  involved  measuring  small  changes  in  potential 
and  distance  as  a  means  of  determining  potential  gradient. 
Tests  were  not  conducted  to  determine  the  accuracy  of  these 
measurements  nor  the  validity  of  the  final  results.   Ref .  12 
states  that  for  a  liquid  electrolyte  it  is  possible  to  obtain 
an  accuracy  of  two  to  five  parts  in  a  thousand  with  respect  to 
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the  theoretical  value  when  using  a  double  probe  of  five  milli- 
meters separation.   The  method  of  determining  the  error  of  such 
an  instrument  was  to  rotate  it  around  a  point  in  the  field  in 
which  case  a  sinusoidal  variation  should  result. 

It  might  seem,  when  first  considered,  that  a  comparison  of 
the  potential  readings  at  corresponding  points  on  adjacent  pro- 
files would  give  an  indication  of  the  overall  accuracy  actually 
obtained.   However,  this  is  not  a  definite  criterion.   In  fact 
it  was  not  uncommon  for  the  measured  potentials  at  correspond- 
ing points  on  adjacent  profiles  to  differ  by  2%   without  produc- 
ing a  noticeable  change  in  results.   The  reason  for  this  is  that 
the  quantities  t    and  ^      which  determine  conditions  in  the  pic- 
ture plane  contain  ratios  of  differences  between  measured  quan- 
tities.  The  measured  values  of  the  various  c  '5  could  there- 
fore differ  from  profile  to  profile  without  necessarily  affect- 
ing the  values  of  F   and  -j- 

A  better  indication  of  the  overall  accuracy  actually  ob- 
tained in  construction  and  measurement  was  the  repeatability  of 
results.   Fig.  33  gives  an  indication  of  the  degree  of  repeat- 
ability of  the  method  and  procedure  as  finally  evolved.   The 
three  layouts  shown  were  distinct  in  every  respect  even  to  the 
use  of  Teledeltos  paper  from  three  separate  rolls. 

The  indication  of  repeatability  does  not,  however,  purport 
to  show  the  accuracy  of  the  method  in  its  degree  of  approxima- 
tion to  the  true  conditions  for  an  infinite  cascade.   The  bus 
bars  were  located  a  finite  distance  from  the  cascade  and  can  on- 
ly approximate  the  conditions  an  infinite  distance  away.   Al- 
though the  error  accepted  with  this  assumption  cannot  be  deter- 
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mined  exactly  with  the  present  equipment,  a  close  approximation 
of  this  error  based  on  a  cascade  of  flat  plates  is  derived  ana- 
lytically in  Appendix  F.   From  the  figure  and  the  data  in  this 
appendix  it  is  evident  that  for  the  geometry  selected,  the  er- 
ror introduced  by  the  location  of  the  bus  bars  is  negligible. 
Only  one  cascade  of  flat  plates  was  investigated,  but  the  ana- 
lytic means  is  at  hand  for  determining  cases  of  different 
solidity  and  stagger  angle. 

The  effect  of  the  foregoing  limitations  on  the  accuracy  of 
K  and   D    cannot  presently  be  determined  directly.   Further 
testing  is  needed  to  provide  statistical  data  for  determination 
of  the  error  as  well  as  means  for  its  reduction. 

The  closely  related  parameter  of  cascade  performance  (9, 
when  compared  to  published  data,  agreed  within  the  limits  of  ±1° 

7.4   CRITIQUE  OF  THE  METHOD  AND  SUGGESTED  IMPROVEMENTS 

The  given  conducting  paper  was  taken  as  fundamental  to  this 
investigation.   The  scope  of  this  work  could  not  include  a 
search  for  nor  an  evaluation  of  various  conducting  mediums. 
However,  removal  of  the  directional  properties  of  the  conducting 
medium  and  reduction  of  the  variation  'in  its  specific  resistance 
should  greatly  improve  the  accuracy  of  the  method.   Perhaps  the 
present  type  of  conducting  paper  could  be  improved.   Other  con- 
ducting mediums  may  be  available.   Considerable  work  has  been 
done  in  France  and  at  Rensselaer  Polytechnic  Institute  in  this 
country  with  electrolytic  tanks.   See  Refs.  12  and  13.   A  sol- 
id conducting  medium  is  believed  to  be  preferable  for  the  meth- 
od herein  described.   The  solid  medium  obviates  such  problems 
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as  leveling  and  polarization.   The  test  set-up  is  much  simpler 
as  only  a  single  template  is  required  for  all  of  the  cascades 
using  that  airfoil. 

If  after  investigation  the  teledeltos  paper  still  proved 
to  be  the  best  solid  medium  available,  allowance  for  the  di- 
rectional properties  could  be  made  by  use  of  a  scale  factor. 
For  the  work  reported  here  a  scale  factor  correction  of  profile 
coordinates  was  not  considered  practical  because  of  the  range 
of  variation  of  specific  resistance  in  a  given  direction.   Sum- 
marizing, even  if  directionality  cannot  be  removed,  variation 
of  specific  resistance  in  a  given  direction  should  be  decreased 

A  minor  suggestion  for  increased  ease  in  cascade  layout  is 
to  change  the  black  color  of  the  backing  side  of  the  sheet  to 
more  clearly  show-up  the  penciled  layout  lines. 

The  layout  table  undoubtedly  would  be  further  modified  if 
this  method  were  to  be  used  commercially.   It  is  considered  de- 
sirable to  bring  the  meridional  ends  of  the  paper  together  to 
establish  an  electrically  continuous  sheet  although  the  origi- 
nal flat  table  produced  satisfactory  results  when  the  approxi- 
mate streamlines  were  used  as  boundaries.   The  cascade  test  at 

(T"  =  1.0  and  p    =  30°  was  done  on  the  original  table.   If  a 
paper  were  used  which  was  sufficiently  long  in  the  direction 
of  the  cascade  axis  the  effect  of  straight-cut  ends  would  prob- 
ably be  small  at  the  center  blade  cutouts.   This  would  be  par- 
ticularly true  for  the  cascades  of  small  stagger  angle. 

If  a  table  be  constructed  of  design  similar  to  the  one 
herein  described,  particular  care  should  be  taken  to  insure 
that  the  table  hinges  be  dimensionally  true.   Their  trueness 
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also  will  require  periodic  checking  because  of  the  dependence  of 
the  electrical  connection  upon  them. 

Perhaps  a  clip-board  type  of  springed  clamp  could  be  design- 
ed for  establishing  electrical  contact  between  bus  bars  and  pa- 
per.  Such  an  arrangement  could  save  time  over  the  bolted  clamp- 
ing system  used  for  these  tests. 

Another  suggestion  for  improvement  is  to  design  the  layout 
table  such  that  the  bus  bars  as  well  as  the  paper  would  be  elec- 
trically continuous.   Although  the  conductivity  of  both  the  cop- 
per bus  bars  and  the  silver  painted  areas  on  the  paper  is  large 
relative  to  the  paper,  their  conductance  would  improve  the  ap- 
proximation of  uniform  conditions  at  large  distances  from  the 
cascade. 

The  method  is  dependent  upon  the  technique  of  the  operator. 
With  an  increasing  number  of  cascades  tested,  the  operator  can 
be  expected  to  improve  in  accuracy  and  increase  in  speed.   The 
angle  of  holding  the  probe,  the  depth  of  probe  penetration,  the 
eyesight  of  the  operator,  and  his  ability  to  position  the  probe 
accurately  are  all  personal  factors  which  will  affect  the  re- 
sults to  a  certain  small  degree. 

The  statistical  curve  fitting  technique  for  determining 
the  value  of  Cj»  probably  could  be  improved.   There  is  a  re- 
quirement for  reduction  of  the  random  error.   The  method  used 
in  this  report  seemed  attractive  because  of  the  similarity  of 
its  equations  to  those  for  perfect  fluid  flow.   It  proved  to  be 
practical  for  solution.   This  same  form  perhaps  could  be  modi- 
fied to  allow  for  the  curvature  of  the  blade  profile  near  the 
trailing  edge. 
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This  method  of  cascade  analysis  has  the  great  advantage  of 
being  much  less  expensive  than  actual  cascade  tests.   The 
teledeltos  paper  is  the  principal  consumable  item.   The  templates 
can  be  manufactured  without  difficulty  or  excessive  expense. 

The  time  saving  of  this  method  over  other  methods  of  cas- 
cade analysis  is  considered  to  be  one  of  its  principal  advan- 
tages.  With  improvement  in  technique  it  was  found  that  a  typical 
cascade  layout  could  be  completed  and  measurements  taken  for  de- 
termination of  turning  angle  in  slightly  more  than  one  working 
day.   This  of  course  depended  to  some  extent  upon  cascade  geom- 
etry.  Calculations  of  turning  angles  from  the  raw  data  were  ac- 
complished in  three  to  four  hours  and  would  be  much  faster  on  new- 
er digital  computers.   The  use  of  the  computer  is  not  required 
for  this  method.   It  was  used  here  principally  to  save  time  and 
thus  permit  the  testing  of  a  greater  number  of  cascades.   The 
measurements  required  for  the  determination  of  velocities  at  60 
stations  on  the  profile  required  an  additional  two  hours.   Cal- 
culating time  to  find  these  velocity  ratios  required  approximate- 
ly  20  minutes  for  each  angle  of  mean  velocity. 

Analyzing  the  method  from  another  viewpoint,  it  experimen- 
tally accomplishes  conformal  transformation  of  an  arbitrary  cas- 
cade onto  a  circle.   When  one  considers  the  well  known  transfor- 
mation which  maps  a  circle  into  a  cascade  of  flat  plates,  he 
might  in  one  sense  consider  that  this  method  determines  an  e- 
quivalent  flat  plate  cascade  for  any  arbitrary  cascade.   The 
performance  of  a  cascade  of  flat  plates  is  defined  theoretically 
by  two  parameters  K  and  &L       .   These  quantities  are  determined 
in  this  method  simply  by  making  the  measurements  required  to 
find  V      and   £-  . 
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This  method  has  greater  capability,  however,  than  simply 
the  determination  of  these  two  parameters.   A  point  to  point 
correspondence  between  blade  profile  and  circle  is  obtained 
which  permits  determination  of  the  local  velocity  ratio  at  any 
point  on  the  blade  profile. 

In  review,  a  method  of  cascade  analysis  has  been  developed 
It  was  shown  to  be  capable  of  predicting,  with  reasonable  accur- 
acy, the  turning  angles  of  perfect  fluid  flow  across  an  arbi- 
trary cascade.   Because  of  the  nature  of  the  present  conducting 
medium  the  accuracy  seems  somewhat  poorer  at  large  stagger  an- 
gles but  no  limits  of  solidity  or  camber  are  known  to  exist. 
The  time  required  to  accomplish  the  procedure  is  short  enough 
to  make  the  method  suitable  for  engineering  application.   The 
method  is  not  expensive. 

Extension  of  the  method  to  the  analysis  of  radial  and  con- 
ical  (mixed  flow)  cascades  is  analytical  once  the  basic  trans- 
formation of  a  rectilinear  cascade  has  been  achieved. 

The  method  could  be  analytically  supplemented  to  permit 
consideration  of  viscous  conditions.   Some  estimate  of  com- 
pressibility effects  could  be  made  using  the  Prandtl-Glauert 
method . 

The  method  of  experimental  conformal  transformation  used 
herein  should  be  a  valuable  tool  when  applied  to  the  many  other 
engineering  problems  which  are  characterized  by  LaPlace's  equa- 
tion.  Examples  are  problems  in  the  theory  of  elasticity,  in 
heat  transfer,  and  in  electric  field  theory. 
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8.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  principal  conclusions  of  this  investigation  are: 

The  method  herein  described  determines  the  ideal  fluid 
flow,  including  turning  angles  and  pressure  distribution 
along  the  blades,  for  any  two  dimensional  cascade.   A 
practical  technique,  partly  experimental  and  partly  ana- 
lytical, of  conformally  transforming  an  arbitrary  cascade 
to  a  corresponding  circle  is  provided  as  the  basis  of  the 
method . 

At  its  present  state  of  development  the  method  appears 
to  determine  turning  angles  to  an  overall  accuracy  of  ±1    , 
but  further  refinement  can  undoubtedly  reduce  this  error. 

The  cascade  analysis  is  accomplished  in  approximately 
one  working  day  and  at  nominal  expense. 

No  limitations  in  the  method  due  to  cascade  geometry 
are  known  to  exist. 
For  further  refinement  of  the  method  it  is  recommend  that: 

A  more  nearly  homogeneous  and  isotropic  conducting  med- 
ium be  used. 

Cascade  layout  and  measurement  procedures  be  mechanized 
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TABLE  1 
INSTRUMENTS 

1.  Potential  supply  -  Heathkit  Battery  Eliminator,  Model  BE-4 . 

2.  Null  meter       -  General  Electric  DC  Micro  Ammeter,  Type 

DO-71,  Model  A29AA1,  coil  resistance  of 
1700  ohms,  detecting  sensitivity  power 
on  the  order  of  8  micro-micro  watts. 

3.  Voltage  divider  -  Electro-Measurements  Inc.  Dekavider,  Mod- 

el RV622,  linearity  of  ±0.0025  percent, 
input  resistance  of  10,000  ohms,  Serial 
13,049. 

4.  Potentiometer    -  Leeds  and  Northrup  Portable  Precision 

Potentiometer,  Number  8662,  Serial 
1318768. 

5.  Voltmeter        -  Weston  DC  Voltmeter,  Model  45,  Serial 

49400. 
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TABLE  II 
AIRFOIL  ORDI NATES  FROM  NACA  ARR.  NO.  L4K22b 
(Stations  and  ordinates  in  percent  of  airfoil  chord) 


Upper 

surface 

station 

Ordinate 

0 

0 

.23 

.97 

.44 

1.24 

.86 

1.66 

2.03 

2.63 

4.43 

4.17 

6.87 

5.44 

9.34 

6.56 

14.34 

8.46 

19.38 

9.98 

24.44 

11.24 

29.55 

12.20 

34.64 

12.90 

39.72 

13.39 

44.82 

13.64 

49.96 

13.65 

55.07 

13.37 

60.17 

12.78 

65.26 

11.92 

70.30 

10.81 

75.32 

9.53 

80.29 

&.04 

85.23 

6.40 

90.17 

4.56 

95.10 

2.57 

100 

0 

Lower 

surface 

station 

Ordinate 

0 

0 

.77 

-.43 

1.06 

-.44 

1.62 

-.46 

2.97 

-.37 

5.57 

-.13 

8.13 

.18 

10.66 

.48 

15.66 

1.08 

20.62 

1.64 

25.56 

2.16 

30.45 

2.60 

35.36 

2.96 

40.28 

3.29 

45.18 

3.58 

50.04 

3.91 

54.93 

4.23 

59.83 

4.48 

64.74 

4.68 

69.70 

4.77 

74.68 

4.67 

79.71 

4.39 

84.77 

3.94 

89.83 

3.10 

94.90 

1.89 

100 

0 
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TABLE  III 
ORDI NATES  FOR  NACA  65-010  BASIC  THICKNESS  FORMS 
(Stations  and  ordinates  in  percent  of  chord) 

Station,  x  Ordinates,  ±y 

65(216)-010  airfoil  combined 
with  y  =  0.0015x 

0  0 

.5  .752 

.75  .890 

1.25  1.124 

2.5  1.571 

5.0  2.222 

7.5  2.709 

10  3.111 

15  3.746 

20  4.218 

25  4.570 

30  4.824 

35  4.982 

40  5.057 

45  5.029 

50  4.870 

55  4.570 

60  4.151 

65  3.627 

70  3.03b 

75  2.451 

80  1.847 

85  1.251 

90  .749 

95  .354 

100  .150 

L.E.  radius  .666 
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TABLE   IV 


ID1  NATES    OF   6"    and    12"   NACA  ARR  NO.    L4K22b  AIRFOILS 


6"  Chord 

(Stations  and 

ordinates  in  inches) 

Upper  surface 

Lower  surface 

Station,  x 

Ordinate,  y 

Station,  x 

Ordinate, 

.0138 

.0582 

.0462 

-.0258 

.0264 

.0744 

.0636 

-.0264 

.0528 

.0996 

.0972 

-.0276 

.1218 

.1578 

.1782 

-.0222 

.2658 

.2502 

.3342 

-.0078 

.4122 

.3264 

.4878 

+.0108 

.5604 

.3936 

.6396 

.0288 

.8604 

.5076 

.9396 

.0648 

1.1628 

.5988 

1.2372 

.0984 

1.4664 

.6744 

1.5336 

.1296 

1.7730 

.7320 

1.8270 

.1560 

2.0784 

.7740 

2.1216 

.1776 

2.3832 

.8034 

2.4168 

.1974 

2.6892 

.8184- 

2.7108 

.2148 

2.9976 

.8190 

3.0024 

.2346 

3.3042 

.8022 

3.2958 

.2538 

3.6102 

.7668 

3.5898 

.2688 

3.9156 

.7152 

3.8844 

.2808 

4.2180 

.6486 

4.1820 

.2862 

4.5192 

.5718 

4.4808 

.2802 

4.8174 

.4824 

4.7826 

.2634 

5.1138 

.3840 

5. 0862 

.2364 

5.4102 

.2736 

5.3898 

.1860 

5.7060 

.1542 

5.6940 

.1134 

6.0000 

0 

6.0000 

0 
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TABLE  IV    (Cont. ) 


12"  Cho. 

rd  (Stations  and 

ordinates  in  inc 

:hes  ) 

Upper 

surface 

Lower 

surface 

Station, 

x  Ordinate,  y 

Station, 

x   Ordinate,  y 

0 

0 

0 

0 

.0276 

.1164 

.0924 

-.0516 

.0528 

.1488 

.1272 

-.0528 

.1056 

.1992 

.1944 

-.0552 

.2436 

.3156 

.3564 

- . 0444 

.5316 

.5004 

.6684 

-.0156 

.8244 

.6528 

.9756 

+.0216 

1.1208 

.7872 

1.2792 

.0576 

1.7208 

1.0152 

1.8792 

.1296 

2.3256 

1.1976 

2.4744 

.1968 

2.9328 

1.3488 

3.0672 

.2592 

3.5460 

1.4640 

3.6540 

.3120 

4.1568 

1.5480 

4.2432 

.3552 

4.7664 

1.6068 

4.8336 

.3948 

5.3784 

1.6368 

5.4216 

.4296 

5.9952 

1.6380 

6.0048 

.4692 

6.6084 

1.6044 

6.5916 

.5076 

7.2204 

1.5336 

7.1796 

.5376 

7.8312 

1.4304 

7.7688 

.5616 

8.4360 

1.2972 

8.3640 

.5724 

9.0384 

1 . 1436 

8.9616 

.5604 

9.6348 

.9648 

9.5652 

.5268 

10.2276 

.7680 

10.1724 

.4728 

10.8204 

.5472 

10.7796 

.3720 

11.4120 

.3084 

11.3880 

.2268 

12.0000 

0 

12.0000 

0 
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TABLE 

V 

COORDINATES  ( 

DF  6"  and  12"  I 

sIACA 

65-(15)10 

AIRFOILS 

6"  Chord 

(Stations  and 

ord: 

Lnates  in  : 

Lncl 

ies  ) 

Upper  surface 

Lower 

sur 

face 

Station,  x 

Ordinate,  y 

Station, 

X 

Ordinate, 

0 

0 

0 

0 

.0059 

.0606 

.0541 

-.0156 

.0181 

.0776 

.0719 

-.0146 

.0438 

.1079 

.1062 

-.0116 

.1122 

.1701 

.1878 

-.0027 

.2558 

.2680 

.3442 

+.0164 

.4033 

.3465 

.4967 

.0351 

.5526 

.4132 

.6474 

.0521 

.8544 

.5229 

.9456 

.0828 

1.1587 

.6079 

1.2413 

.1085 

1.4643 

.6746 

1.5356 

.1309 

1.7709 

.7254 

1.8291 

.1494 

2.0780 

.7603 

2.1220 

.1667 

2.3853 

.7850 

2.4147 

.1789 

2.6934 

.7944 

2.7672 

.1911 

3.0000 

.7886 

3.0000 

.2042 

3.3066 

.7669 

3.2934 

.2186 

3.6120 

.7307 

3.5880 

.2332 

3.9160 

.6796 

3.8640 

.2474 

4.2184 

.6188 

4.1816 

.2560 

4.5191 

.5486 

4.4809 

.2569 

4.8181 

.4675 

4.7819 

.2489 

5.1152 

.3763 

5.0848 

.2294 

5.4114 

.2761 

5.3886 

.1892 

5.7070 

.1622 

5.6930 

.1222 

6.0000 

.0090 

6.0000 

-.0090 

LEADING  EDGE  RADIUS  =  .03996  inches 
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TABLE  V  (Cont. ) 


12"  Chord 
Upper  s 
Station, , 

0 

.0116 

.0362 

.0876 

.2245 

.5116 

.8066 

1.1053 

1.7088 

2.3174 

2.9287 

3.5417 

4.1559 

4.7707 

5.3868 

6.0000 

6.6131 

7.2241 

7.8321 

8.4367 

9.0382 

9.6362 

10.2304 

10.8228 

11.4141 

12.0000 


(Stations  and 
urf ace 

x  Ordinate,  y 

0 

.1213 

.1553 

.2159 

.3401 

.5359 

.6930 

.8264 

1.0459 

1.2158 

1.3492 

1.4507 

1.52C6 

1.5700 

1.5888 

1.5771 

1.5337 

1.4614 

1.3592 

1.2375 

1.0971 

.9351 

.7527 

.5522 

.3245 

.0180 


ordinates  in  inch 
Lower 
Station, . 

C 

.1082 

.1438 

.2124 

.3755 

.6864 

.9934 

1.2947 

1.8911 

2.4826 

3.0713 

3.6583 

4.2440 

4.8293 

5.5344 

6.0000 

6.5868 

7.1759 

7.7679 

8.3633 

8.9618 

9.5638 

10.1696 

10.7772 

11.3859 

12.0000 


es  ) 

surface 
x  Ordinate,  y 
0 
-.0313 
-.0293 
-.0233 
-.0053 
+.0328 

.0702 

.1042 

.1655 

.2172 

.2618 

.2989 

.3334 

.3578 

.3822 

.4083 

.4373 

.4664 

.4949 

.5121 

.5139 

.4977 

.4587 

.3784 

.2443 
-.0180 


LEADING  EDGE  RADIUS  =  .07992  inches 
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TAB'LE  VI 
TEST  PROGRAM 


A.   Turning  Angles: 
1.   First  series 


Profile  -  Mutterperl  (NACA  ARR  No.  L4K22b ) 
Tests      (a)  £>  =  45°;  (T  =  1.0 
(b)  (3  =  450;  (T  =  1.5 


2.   Second  series: 

Profile  -  NACA  65-( 15)10 
Tests  -  Group  one  (six  inch 


1.5  (twleve  inch  chord) 

1.5  (six  inch  chord) 

1.5  (twelve  inch  chord) 

1.5  (twelve  inch  chord) 


B.   Velocity  Distribution: 

Profile  -  NACA  65- ( 15)10  (twelve  inch  chord) 
Test  -  p   =  35°;   (T  =  1.5 
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TABLE  VII 
SAMPLE  DATA  FOR  TURNING  ANGLE  DETERMINATION 


Upp 

er 

Sur 

face 

Station 

6 

(Percent ) 

3b 

71.88 

c 

71.88 

d 

71.89 

4 

71.90 

4a 

71.89 

b 

71.88 

46 

34.91 

a 

34.70 

b 

34.47 

c 

34.26 

d 

34.00 

47 

33.77 

a 

33.52 

b 

33.28 

c 

33.01 

d 

32.76 

48 

32.50 

a 

32.20 

b 

31.93 

c 

31.59 

d 

31.24 

49 

30.68 

Center  Tail 

30.52 

Lower 

Sur 

face 

Station 

£ 

(Percent,) 

46 

30.69 

a 

30.61 

b 

30.55 

c 

30.49 

d 

30.42 

47 

30.37 

a 

30.29 

b 

30.24 

c 

30.18 

d 

30.12 

48 

30.08 

a 

30.03 

b 

30.01 

c 

29.99 

d 

30.00 

49 

30.18 

ITx/eo 

• 

Station. 

.       £ 

(Percent) 

0 

98.53 

1 

95.74 

2 

92.97 

3 

90.24 

4 

87.58 

5 

84.84 

66 


TABLE  VIII 
SAMPLE  DATA  AND  INPUT  VALUES  FOR  DETERMINATION  OF  VELOCITY  RATIOS 


Upper  Surface 


Lower  Surface 


Percent 

<5 

a£ 

,&  , 

,   6     , 

±E 

£T 

Chord 

(Percent ) 

(Volts) 

(Percent) 

(Percent ) 

(Volts) 

(Percent ) 

From  L.E. 

0 

68.38 

67.84 

14.17 

2.06 

71.31 

10.22 

2.353 

65.30 

34.86 

7.963 

4.17 

71.75 

3.58 

0.824 

63.67 

29.33 

6.752 

6.25 

71.85 

- 

- 

62.22 

27.45 

6.319 

6.33 

71.75 

2.57 

0.592 

61.00 

22.97 

5.288 

10.42 

71.53   ' 

4.77 

1.098 

59.85 

21.70 

4.995 

12.50 

71.23 

6.06 

1.395 

58.74 

19.71 

4.537 

16.67 

70.40 

8.91 

2.051 

56.71 

19.29 

4.441 

20.83 

69.29 

10.69 

2.461 

54.65 

18.75 

4.316 

25.00 

67.99 

13.14 

3.025 

.52.81 

18.70 

4.305 

29.17 

66.51 

13.93 

3.207 

51.00 

16.30 

3.752 

33.33 

64.88 

15.76 

3.628 

49.26 

16.23 

3.736 

37.50 

63.06 

18.23 

4.197 

47.56 

15.75 

3.626 

41.67 

61.21 

18.43 

4.243 

45.95 

16.25 

3.741 

45.83 

59.18 

18.73 

4.312 

44.34 

15.77 

3.630 

50.00 

57.15 

20.31 

4.675 

42.74 

15.77 

3.630 

54.17 

55.06 

20.12 

4.632 

41.22 

14.70 

3.384 

58.33 

52.98 

19.64 

4.521 

39.81 

13.15 

3.027 

62.50 

50.87 

19.53 

4.496 

38.44 

13.08 

3.011 

66.67 

48.81 

18.30 

4.213 

37.14 

12.17 

2.802 

70.83 

46.68 

18.98 

4.369 

35.94 

11.25 

2.590 

75.00 

44.57 

19.40 

4.466 

34.79 

11.54 

2.657 

79.17 

42.52 

18.45 

4.247 

33.72 

10.00 

2.302 

83.33 

40.39 

18.34 

4.222 

32.77 

9.22 

2.122 

87.50 

38.23 

19.31 

4.445 

31.89 

8.21 

1.890 

91.67 

36.10 

20.55 

4.731 

31.05 

6.80 

1.565 

93.75 

34.91 

20.02 

4.609 

30.69 

5.94 

1.367 

95.83 

33.97 

19.85 

4.570 

30.37 

5.35 

1.232 

97.92 

32.50 

20.70 

4.765 

30.08 

3.71 

0.854 

AZ  =  .0405  inches 
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TABLE   IX 
INPUT   VALUES   FOR  DETERMINATION  OF   TURNING  ANGLES 


<r 

Pearees 

O/                  oz                   Cf 

Percent             Percent             Percent 
MUTTERPERL  AIRFOIL 

~TX 

Percent 

5 

l,nche§ 

1.0 
1.5 

45 
45 

61.86                 39.42                 40.29 
62.61                 38.02                 38.54 

2.870 
2.825 

6.00 
4.00 

NACA   65- (15) 10 


1.0 

20 

61.65 

40.54 

40.63 

1.0 

25 

61.14    . 

40.32 

40.50 

1.0 

30 

62.14 

39.01 

39.30 

1.0 

35 

62.34 

40.82 

41.26 

1.5 

20 

71.38 

30.26 

30.36 

1.5 

25 

61.63 

39.70 

39.83 

1.5 

30 

72.26 

30.03 

30.44 

1.5 

35 

71.90 

29.99 

30.49 

3.060 

6.00 

3.072 

6.00 

3.259 

6.00 

2.882 

6.00 

2.943 

8.00 

3.019 

4.00 

2.761 

8.00 

2.761 

8.00 
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TABLE  X 
SENSITIVITY  OF  RESULTS  TO  VARIATION  OF  INPUT  VALUES 


Quantity 

Value 

AQuantity 

^Quantity. 

Quantity 

A® 
atfl(5n=   35° 

A<9 
atofm=   60° 

C, 

62.91 

+6 

.001 

+    .012° 

+.026° 

&. 

38.58 

+4 

.001 

+1.105° 

+.641° 

Of 

39.68 

+4 

.001 

-1.104° 

-.647° 

30.80 

+  3 

.001 

+    .029° 

+.016° 
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Figure  4  -  Transformations  Between  The  Real  and  Picture  Planes 
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Figure  5  -Simple  Two  Dimensional  Flow3 
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Figure  6  -Correspondence  Between  The  Real  and  Picture  Planes 
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Figure  20o-  Cutting  knife  modifications, 
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APPENDIX  A 
THE  ANALOGY 


in. 


In  potential,  two  dimensional,  ideal  fluid  flow  there  is  a 
complex  velocity  potential  or  =  <p   +  Ity      such  that  the  Cauchy- 
Riemann  differential  equations  are  satisfied 

and 

Then  the  complex  velocity  may  be  written 


I 


Al 


d^    tt  ^r  i^r  ^f 


A2 


Considering  the  flow  of  a  steady  current  in  a  flat  sheet  of  uni- 
form resistance  per  unit  square  y    . 


••  % 


A3 


where  <?n   is  the  current  line  intensity.   Let  1  be  a  total  cur- 
rent function  analagous  on  the  sheet  to  the  stream  function  in 
the  ideal  fluid  flow  such  that 


i»    =-~  and  I     9*L     . 


A4 


It  follows: 


-p^L  =  _^iL  and  v^-=^- 


A5 


which  indicates  that  E  and  I  also  satisfy  the  Cauchy-Riemann  dif- 
ferential equations.   Equations  A5  may  be  used  to  obtain 


E+lvI  =  f  (***</)  =f(2)   , 


A6 
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or  if  we  divide  by  the  maximum  potential 

E*=Ki*  ,  A7 


f*+^  =  £^  =*,(*)  . 


Hence  we  may  list  the  following  analogies: 

2   Plane  o*  Plane 


A8 


A9 
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APPENDIX  B 
DERIVATION  OF  EQUATIONS  FOR  CONFORMAL  TRANSFORMATION 

This  appendix  develops  the  equations  needed  to  achieve  con- 
formal  transformation  of  the  cascade  of  blades  in  the  z  plane  to 
the  circle  in  the  *f  plane.   It  also  defines  the  quantity  -£p£- 
which  is  important  in  later  work  for  relating  velocities  in  the 
picture  plane  to  those  in  the  real  plane. 

Conformal  transformation  is  possible  between  the  real,  z 
plane  and  the  picture,  5    plane  if  there  exists  a  function  in 
the  picture  plane  of  the  real  plane  complex  variable  z  such 
that  for  every  assigned  value  of  the  variable  the  function  has  a 
definite  value  or  system  of  values  and  also  a  definite  derivative. 
The  advantage  of  this  definition  is  that  it  does  not  require  the 
existence  of  an  analytical  expression  of  the  function. 

It  was  shown  in  the  main  text  that  any  point  on  the  real 

plane  may  be  designated  by  its  coordinates  x,  y  or  by  the  values 

of  O  and  ^*r  which  characterize  it.   See  also  Appendix  A.   Fur- 

thur,  the  parameter  X  —       i  $i\      was  shown  to  be  the  unique  and 

•5  v"3x  )<x> 
invariant  quantity  which  must  be  maintained  to  achieve  conformal 

transformation. 

The  basic  potential  relations  for  a  point  of  strength  m  are 

p  =^rlnK      and   ^  =  m.  Q  B1 

where  m  is  real..   The  point  is  called  a  source  if  m  is  positive 
and  a  sink  if  m  is  negative.   R  is  the  radial  distance  from  the 
point  and  Q     is  the  angular  coordinate  with  respect  to  the 
point  as  origin. 
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In  the  picture  plane  of  Fig.  4,  let  sources  of  strength  m 
be  located  on  the  real  axis  at  5  ="0,-0,^  and  let  sinks  of 
equal  strength  be  located  on  the  real  axis  at£=b»Q/k  . 
Then  the  complex  potential  in  the  picture  plane  is 


or  = 


=  _ni[n_bj-l_  +  j2L[n  t  +  ty   m  P+ilp    • 


B2 


zn  b— J5     T^"'jf-% 

Now  any  point  in  the  picture  plane  can  be  designated  by  its  coord- 
inates J  > 1  or  by  the  values  of  f  and  y  which  characterize  it. 
Combining   the  factors   of    eq.    B2  , 


B3 


On  the  circle 


=  a  , 


or  -  m  (n  Qe*  [Qe**  +  Qe"j9  +  b//(X  +  %- 


B4 


and  with    tyfa  =  e 


-  ** 


ttr  -  m  /  n  cosh  y  +  g  cos  6 
in     cosh  T  -  a  c  os  0 


B5 


W 


e  may   let  a   =   1,    giving 


or  -  rn  \r\  ^os^y  +  cos  6 
~"  zrr       coshr  -  cos© 


s  p  +  2^ 


B6 


which  shows  that  the  circle  is  a  line  of  constant  W 
If  1  and  2  are  the  points  indicated  in  Fig.  4 


B7 


V  _  (D  -  m  In  cos/ifr*  1 

By  Appendix  A,  O  corresponds  to  C   ;  so  we  may  impose  the  rela- 


tion 


f  =  /c,-k26    , 


B8 
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where     k,     and     k,    are   arbitrary  constants.      It   follows   that 


fc-f/  =  M£f-<54) 


and   if 
then 


m-ak^-g) 


CO 


TO.        J I     -     «p»     —     K*  \Ci      CaJ 


B9 

BIO 

Bll 


This    establishes    the   existance   of   a   function  relating   the   z   and 
5    planes   conformally. 


Eq.    B7   may   be    solved    to   give 


tf  =  cosh 


-/ 


•m 


7%^ 


+  L 


0**  - 


B12 


=     cosh 

e   ™    -  / 

Knowing  '2f  ,  the  locations  of  the  sources  and  sinks  in  the  pic- 


1TV 


ture  plane  are  determined.  For  any  point  P  on  the  circle  and  a 
corresponding  point  P'  on  a  blade  profile,  let  there  be  a  quan- 
tity £  such  that 


B13 


By  substitution  of  eqs .  B6  and  B7  into  eq.  B13  we  obtain 


cosh  V  +  cos© 


I     cosn  jrn>jc7       t    coshZ+l 
T      Z  i„  cosh?  +  I 


B14 


(r? 


coshtf-  / 


and   solving   for    cos  0  : 


cos  5  =  cosh  t 


[coshV-lJ        ~~    I 


cosh*-)/        +    / 


=  coshT 


e  -  I 


etrr^-/)  + 


B15 


The  derivative 


dl 


d2 


is  the  operator  which  relates  infin- 
itesimal vectors  in  the  picture  plane  to  corresponding  vectors 
in  the  real  plane,  and  vice  versa.   In  particular, 
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_  c/ur 
~di     ~   d$ 


V,  =* 


c(5 
d2 

=  v$ 

d* 
d2 

B16 


but  recalling  our  parameter  -r   this  may  be  written  t 


Vf  =  Vj 


d2 


if 

d5 


B17 


Nov  V5  is  known  in  the  picture  plane  and  is  given  by  eq.  C5  of 
Appendix  C. 

The  numerator  of  eq .  B17  can  be  written  as 


df 


di 


&i 


e,-£»    • 


The  denominator  of  eq.  B17  can  be  written  as 


B18 


dS 


Sir\9  coshtf 


(cosh^-co5z©)ln-^lfn 


cosher 


B19 
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APFENDIX  C 
DERIVATION  OF  FLOW  EQUATIONS 

This  appendix  derives  the  equations  required  to  develop  an 
arbitrary  flow  about  the  circle  pp  a  . 

The  sources  and  sinks  used  to  accomplish  transformation  can 
be  considered  to  create  an  axial  flow  velocity  \/a  such  that 


)T1  =  VaS  .   On  the  circle  5 
locity.   From  eq.  B3 


=  a   this  flow  produces  a  ve- 


1     ^ P  .  Vo  -      2^5      coshtsinO 


R  7© 


77     '  coshZV-cosze 


Cl 


The  basic  velocity  and  stream  potential  relations  for  a  point 
vortex  of  strength  J\,   are 


V-^-9 


and  y^-J^/ntf, 


C2 


where  7\.  is  positive  if  counterclockwise. 

If  vortices  of  strength  J\   and  with  direction  indicated  in 
Fig.  7  be  placed  in  the  picture  plane  in  addition  to  the  sources 
and  sinks  the  resulting  complex  potential  is 


C3 


and  the  corresponding  velocity  on  the  circle  is 
If  \/a=f-  and  V*  = 


7T    "  C0&h21f-C0S29 

.A 


7?   COSh^^-COS^e 


C4 


C5 


Since  the  value  of  m  is  determined  in  the  process  of  transforma- 
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tion,  from  m  =  5  KJ — ^-]   >  "the  arbitrary  choice  of  a  mean  veloc- 
ity direction  oCm        determines  the  strength  -A.   .   Then  if  ad- 


ditional vortices  of  strength 


r. 


and  with  direction   indicated 


in  Fig.    7   be   placed   in   the  picture  plane,    the   resulting   velocity 


on   the   circle 


i 


=  a     is : 


\fa  n   V  cosher- cosze  )   zn\coshZi-cosZ&)>  uo 

At  the  point  on  the  circle  where  6  =  Of ,  \jq  =•  0   •   Solving  for 


JL 


r   =  tan  oTm    C0SL°r  -  -M&l    . 


Zm  '  7%s  ~  vw"  ^m    coshr       sinhZ 


C7 


The  circulation  given  by  this  equation  insures  satisfaction  of  the 

Kutta  condition. 

The  velocity  triangle  is  now  complete. 


ton  qC<  =  tcmoCm  +  -^f- 


ton  oCz  -  tan  of™  -  -£- 


tan(dC-&)  = 


r 


ce 


C9 


cio 


=  tan  8. 

Of  aerodynamic  interest  is  the  ideal  angle  of  attack.  At  this 
angle  the  front  stagnation  point  is  at  the  leading  edge  of  the 
airfoil.   If\^=0  IccS^Ol   equation  C6  may  be  solved  to  give, 


tnn  rSi  -    sin  Or  -sinQi,  _ 

For  the  determination  of  lift  coefficients, 


cn 


C12 


From  eqs.  C8  and  C9  ,    Ts  Va^(tQJl  cC,  —  tQDoCi ) 


C13 
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Substituting   and   solving, 

,      =    zVrnVaS  fmoCi  _  tar)oCl)  C14 

Hence , 


CIS 
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APPENDIX  D 
DEFINITION  AND  EXPLANATION  OF  THE  TERM  "SPECIFIC  RESISTANCE" 


Consider  an  element  of  the  conducting  medium  of  length  A   1 
and  width   A  W.   Assume  constant  thickness  h  and  resistivity  p  . 
Let  ZW  be  the  potential  drop  along  the  length  Al  of  the  ele- 
ment, and  Ai  be  the  corresponding  current. 


&l 


Ohm's  Lav:  may  be  written  in  the  form: 

AV^L  =  p  Ai/AA   =(/%)(^w)  where  A/^  =  K^W). 
The  constant  of  proportionality  P/h   is  defined  as  the 
specific  resistance  V  ,  of  the  conducting  medium.   Solving  the 
above  equation  for  this  term: 

■v  =  (PA)  =  &/*i  (aWai)  rffi8  in„-  =  ^ • 

4 

The  net  resistance  of  the  element  may  be  written: 

r  =  aV/a?  =  vf%,). 

Thus  if  the  element  is  square  (Al  =AW),  its  net  resist- 
ance R,  will  be  numerically  equal  to  the  specific  resistance  *V  , 
of  the  medium,  regardless  of  the  size  of  the  square. 
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APPENDIX  E 
DETERMINATION  OF  £> 


Let  UTa  =AzPa  =  £+i&  El 


and  bottom  of  a  wedge,  Fig.  El,  where  pa  —   23 —  • 


be  the  complex  potential  for  conditions  symmetrical  on  the  top 

/ 

1 

2-rr 

The  sides  of  the  wedge  are  two  separate  lines  which  join  at  the 
wedge  vertex  to  form  a  single  line  bisecting  the  obtuse  angle 
formed  by  the  vertex.   Along  each  line  ^  =  0;  so, 

6a  =  n  T  E3 

Similarly  let       W~b   =  6  l?b  E4 

be  the  complex  potential  for  conditions  where  a  single  line  of 
constant  \£   coincides  with  the  top  and  bottom  boundaries  of  the 

wedge.   Here   Pk-J^    =  &   Pa  ~  P  E5 

and  as  above,  £,    =  ^Of  E6 

on  the  vyedge  profile.   The  +  sign  is'  for  the  upper  side  of  the 
wedge  and  the  -  sign  is  for  the  lower  side  of  the  wedge. 

Superposition  of  these  two  complex  potentials  results  in 

Let  the  wedge  approximate  the  trailing  edge  of  the  airfoil 
being  investigated.   Let  r  be  the  straight  line  distance  of  a 
point  on  the  profile  from  the  trailing  edge.   With  the  trailing 
edge  at  the  point  r  =  0,  eq.  E7  will  result  in  O    being  zero  at 
r  =  0;  so  with  the  addition  of  a  constant,  the  resultant  equation 

is,  £d  =  CT  +Ari?  ±  drp  .  E8 
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Desired  is  the  determination  of  the  constants  Cjji  ,   A,    B, 
and  p,  such  that  the  curve  of  the  above  expression  will  as  close- 
ly as  possible  approximate  the  measured  values  of  £/    in  the 
neighborhood  of  the  trailing  edge  of  the  airfoil  and  permit  de- 
termination of  o«  .   The  constant  p  is  determined  by  a  plotting 
procedure  described  later  in  the  appendix.   The  determination  of 
0  y     is  accomplished  by  the  method  of  least  squares  which  is 
mathematically  expressed  by  making  the  sum  of  the  squares  of 
the  deviations  of  the  experimental  points  from  the  analytical 
curve  a  minimum: 

A  =  Z-j\&i  ~  MjJ  =  minimum  .  E9 


This  requires  that , 


14-  n  -*A  -1A 


o*=«„-Tiriitfp,      avau.fi&rr  Mz^t^Si 

c*i  l~\  i=i 


E1C 


The  resulting  set  of  equations  is 

Q(lUr  +  <XiZf\   +  a,3  8  =  M, 

where 


E12 


and  N  is  the  number  of  experimental  points. 

The  matrix  solution  for  C-j*    is  obtained,  yielding 

r  ^a2Q.,8Qi5  +  a^a35  -  azi  -a,2a33  -  a,3  aiZ 
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For  this  investigation  a  sketch  of  the  trailing  edge  of  the 
airfoil  was  drawn  to  a  scale  ten  times  that  of  the  6  inch  chord 
templates  used.   Straight  line  distances  were  measured  on  this 
sketch  from  the  trailing  edge  of  the  airfoil  to  the  locations  a- 
long  the  top  and  bottom  surface  where  potential  readings  were 
taken.  A  curve  was  plotted  of  Cu   and  o„   versus  these  distances 
r;  See  Fig.  E2.  From  this  curve  the  values  of  Ou  and  Og     were  ob- 
tained at  equal  values  of  r.   On  log-log  coordinates  (Ou  "~  w  ) 
was  plotted  as  a  function  of  r.   The  slope  of  this  straight  line 
is  the  coefficient  p.   See  Fig.  E3. 
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APPENDIX  F 
APPROXIMATION  OF  THE  ERROR  DUE  TO  THE  BUS  BAR  LOCATION 

In  adapting  equipment  to  provide  the  electrical  analogy  to 
flow  through  cascades,  it  was  necessary  to  locate  the  bus  bars 
a  finite  distance  from  the  cascade.   Thus  the  location  of  the 
bus  bars  introduced  an  error  in  the  measured  data.   The  deter- 
mination of  this  error  cannot  readily  be  made  when  the  airfoils 
of  the  cascade  have  camber  and  thickness.   However,  a  good  ap- 
proximation of  the  error  may  be  found  by  considering  a  cascade 
of  flat  plates. 

Tests  of  this  investigation  were  conducted  under  conditions 
equivalent  to  straight  through  flow  with  no  circulation.   The 
same  conditions  may  be  achieved  over  a  flat  plate  cascade  by 
imposing  a  potential  flow  composed  of  sources  and  sinks  onto  the 
conformal  map  of  the  cascade. 

THE  CASCADE  OF  FLAT  PLATES 


Considering  y -   f(5)      as  a  complex  potential  function, 


then 


i-sfSP- 


zn 


ln^±   -In 


Fl 


gives  the  mapping  function  between  a  circle  with  the  radius  a 
in  the  $  plane  and  a  staggered  lattice  composed  of  straight 
airfoils  in  the   ^T   plane  (Fig.  Fl). 

The  first  term  of  the  above  equation  can  be  considered  as 
originating  from  sources  and  sinks  of  equal  strength  located  as 
shown  in  the  figure,  while  the  second  term  originates  from 


119 


vortices  of  equal  magnitude  superimposed  on  the  sources  and 
sinks.   The  sense  of  the  vortices  is  shown  by  the  direction  of 
the  arrows  in  Fig.  Fl. 

The  spacing  of  the  airfoils  along  the  cascade  is  a  con- 
stant s  and   p   if  the  stagger  angle  of  the  blades. 

The  parameter   b  is  connected  with  the  chord  spacing  ratio. 

In  order  to  establish  the  relation  between  them,  it  is  necessary 

to  calculate  the  chord.   By  settingd^  =  0  for  5=  5L(L,F;y.  Fl)  and 

a  5 

$-$T(TiFiq.F  I)    which  points  correspond  to  the  leading  and 
trailing  edges  of  the  blade,  it  can  be  shown  that 


j-^-^cosp  In 


\Vb*  +  2/bz cos%p  +  1  )*     +  z/b  cos  ft 


-'/b2 


F2 


+  sinp  tan1 


%  sin  p 


(&♦  +  wcofp+\y 


For  any  chosen  stagger  angle  8  ,  the  solidity  (T=  c/s  may  be 
computed  for  various  assumed  values  of  b.   By  plotting  these  re- 
sults, it  is  possible  to  determine  b  graphically  for  any  choice 
of  (T  •   Table  F  I  gives  a  set  of  such  corresponding  values  for 
a  p    of  45°. 

In  the  present  example  the  solidity  was  arbitrarily  chosen 
as  unity.   The  parameters  of  the  cascade  geometry  are  now  fixed. 

In  order  to  obtain  the  point  to  point  relation  between  the 
'V  and  ^  plane  it  was  necessary  to  reduce  equation  Fl   to  its 
real  and  imaginary  parts.   After  a  simple  but  somewhat  lengthy 
calculation  there  is  obtained  the  form, 
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i=xV^=5f^ln 


{(r  +#)cos  0  +  Zacoshr]  +  2  {(<#  -  y-)sm9} 


((r  +  a>>o5e -Zacoshr}  +  j  {(a>  - r)s*nfi/ 


-t-Jtan$  in 


{(r-Q>)cos9^2a^^^(^^r)5in6J 


j(r-  <ft)cos0-2<n  smhrj  +  a  [(#.  +r)  sm©} 


Where     5  =  Y*6 


?<9 


tr 


a  nd   b  —  C  • 


Fla 


THE  STRAIGHT  THROUGH  FLOW 

< 

The  complex  potential  function  necessary  to  produce  a 
straight  through  flow  with  no  circulation  is 


3-b 


F3 


where  \(^S    represents  the  strength  of  the  sources  and  sinks  in 

the  5  plane.   In  the  y-  plane  V^  is  the  velocity  at  infinity 

and  s  is  the  cascade  spacing.   The  product  VqS   ,  therefore  rep- 

sents  the  quantity  of  flow  passing  between  two  adjacent  blades. 

The  deviation  of  the  velocity  along  any  streamline  from  the 

value   Vet  ^s  a  measure  of  the  error  in  locating  the  bus  bars  a 

finite  distance  from  the  cascade.   If  the  ratio  Vfl  Vfr  along 

Va 
any  streamline,  at  a  position  along  the  streamline  corresponding 

to  the  location  of  the  bus  bars,  is  negligible  then  the  error  due 

to  the  bus  bars'  location  is  negligible. 


121 


The   stagnation   streamline  was   chosen   for   investigation   as 
it   has    the  most   radical   changes   in  velocity   and   direction.      Cal- 
culations   for   this   streamline   in   the    o      plane   occur  where   Q~Oi 
ft   and     b  ^f^  Q     •      Setting   Q=   1   for   convenience,    and  placing 
0  =  TT,    the  mapping   equation   reduces    to: 


Z1T 


(r  +  yc)-Zcosh% 


Zcoshr+(tiVr\ 


Fib 


\r-Vf)-Zsinh7 


ZslnhZ+(r-  W)_ 


In  order  to  facilitate  visualization  of  perpendicular  distances 
from  the  cascade  centerline  a  simple  rotation  of  axes  in  the  'V 
plane   is    introduced, 

^zzx  +  iy      ushere    %  =  X'ccsQ  -tfsinfi 

and      y=  x'sinp  +  ifcosQ     . 


F4 


The  complex  velocity  in  the   "^r   plane  is, 


*v 


With  the  selections  of  0  =  7?,  Q  =  1.00  and  6 
mines  'X'   for  this  case, 


F3a 


=  45  which  deter- 


V*.  I 

""ff-if   +  ita^2T(r4-/) 

Va      cos£ 

(V-/)*  4  taT>h*r(r,'+ir 

F3b 


The  magnitude  of  this  ratio  does  not  change  under  a  rotation 
of  the  axes. 

The  value  of  f  in  the  5   plane  was  varied  from  1.0  to  b. 
This  gave  corresponding  points  in  the  \     or  **y    plane  from 
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equations  Fib  and  F4  respectively.   At  the  same  time  the  ra- 
tio  a   ^  was  computed  through  the  use  of  equation  F3b. 

The  values  of  x/s  and  corresponding  values  of  Y3  Vfr 

Va 
for  the  stagnation  streamline  of  the  cascade  are  presented  in 

Table  FII  with  T    as  the  argument.   Fig.  F3  provides  a  graphic 

display  of  the  results.   The  cascade  in  the  figure  is  drawn 

to  scale. 

The  working  surface  as  finally  evolved  had  bus  bars  located 

16  inches  from  the  cascade  centerline.   Two  sizes  of  profiles, 

6  and  12  inch  chord,  were  used  in  the  investigation.   The 

solidity  for  the  flat  plate  cascade  has  been  selected  as  unity 

thereby  limiting  the  values  of  the  ratio  x/s  to  16/12  and  16/6. 

These  two  cases  are  shown  in  Fig.  F3. 
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TABLE   FI 
CORRESPONDING  VALUES  OF      b     and    (T    at   6  =   45( 


b 

(T 

at  t 

)    =  45° 

2.00 

.6560 

1.95 

.6739 

1.90 

.6930 

1.85 

.7133 

1.80 

.7350 

1.75 

.7586 

1.70 

.7837 

1.65 

.8111 

1.60 

.8413 

1.55 

.8737 

1.50 

.9100 

1.45 

.9506 

t 

1.40 

.9966 

1.35 

1.0491 

1.30 

1.1107 

1.25 

1.1846 

1.20 

1.2765 

1.15 

1.3969 

1.10 

1.5699 

1.05 

1.8716 

1.00 

00 
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VQ 


T 
±.2959 


.9131 


.3535  cascade  edges 


TABLE  FII 
CORRESPONDING  VALUES  OF  -|-  AND   Va J  ^»   at   (T  =  1.0  and  ^ 

r 

1.02 

1.04  +.3061  .8289 

1.06  ±.3155  .7211 

1.08  ±.3274  .6724 

1.10  ±.3394  .5993 

1.12  +.3519  .5324 

1.14  ±.3650  .4793 

1.16  ±.3794  .4113 

1.18  ±.3948  .3537 

1.20  ±.4111  .3186 

1.22  ±.4345  .2636 

1.24  ±.4503  .2227 

1.26  ±.4732  .1846 

1.28  ±.4994  .1505 

1.30  ±.5309  .1189 

1.32  +.5684  .0997 

1.34  ±.6190  .0639 

1.36  ±.6909  .0494 

1.38  •  ±.8220  .0168 

1.39  +.9743  .0055 


=  45 
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